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TABLE  OF  SYMBOLS 

Symbol  Quantity  Units 

A  =        cross  sectional  area  of  metal  ft 

m 

2 

A  =        cross  section  area  of  fluid  ft 

f 

a  =        inner  radius  of  pipe  wall  ft 

a  =        half  width  of  channel  (slab  geometry)        ft 

a  =        b,b2A.     »   used  in  equations  (IV-29) 

and  (IV-31)  hr-1 

Bf  =        ha/Kf  =  Fluid  Biot  number  defined 

by  equation  (V-33) 

Bf*  =        2ha/KF  =  Fluid  Biot  number  defined 

by  equation  (V-62) 

B  =        hb/K       =  pipe  wall  Biot  number 


m 


c 


defined  by  equation  (V-33) 


b  =         outer  radius  of  pipe  ft 

b  =        thickness  of  wall  (slab  geometry)  ft 

b.  =        *"/   /fOfAj.  =  parameter  defined  by 

equation  (IV_2)  hr 

b  =  /    P    CrnArn  =  parameter  defined 

by  equation  (IV_5)  hr 

C£  =         Fluid  specific  heat  /lbm°F 


m  =        pipe  wall  specific  heat  BtuA,       o_ 


C^  =        coefficient  defined  by  equation  (IV-57), 

i  =  0,  1 

C^  =         coefficient  defined  by  equation  (IV-67), 

1  i  =  0,  1 


h* 

h(f,n;t) 

In(x) 

I 
n 


Jn(x) 


Kma 

Kmt 

K 
m 

Kfa 


Lft 


K 


ki 


*i 


Kn(x) 


P 

q 


heat  transfer  coefficient  at  fluid-wall 
interface 

heat  transfer  coefficient  defined  by- 
equation  (II-  15) 

heat  transfer  coefficient  defined  by- 
equation  (II- 9) 

Function  defined  by  equation  (IV-45) 


modified  Bess  el  function  of  first  kind 
of  order  n 

complex  integral  defined  in  Appendix  A 
equation  (5) 

Bessel  function  of  first  kind  of  order  n 

pipe  wall  axial  conductivity 

pipe  wall  transverse  conductivity 

pipe  wall  isotropic  conductivity 

fluid  axial  conductivity 

fluid  transverse  conductivity 

fluid  isotropic  conductivity 

coefficients  defined  by  equation 
(IV-40),   i  =  0,  1,2,3,4 

coefficients  defined  by  equation 
(IV-42),   i  =  1,2,3,4,5 

modified  Bessel  function  of  second 
kind  of  order  n. 

perimeter 

Is/*)"2 

radial  coordinate 


B*Vft2hr°F 
Btu/ft2hr°F 
Btu/ft2hr°F 


Btu/ft  hr°F 
Btu/ft.hr°F 
Btu/ft.  hr°F 
Btu/ft.  hr°F 
Btu/hr  ft  °F 
Btu/hr  ft  °F 


ft 
ft 
ft 


-1 


Laplace  transform  parameter 


hr 


"I 


.11  On 

T                           =        PxPe  wall  temperature  F 

T*  =         radially  averaged  pipe  wall 

temperature  F 

t                             =         Time  hr 

u.                          =        temperature  functions  used  in 

Section  III  F 

V                          =        fluid  velocity  ft/hr 


o 


W(x,  y,  t)  =        temperature  function  defined  by 

equation  (111-21)  WF 

W(r,x,  t)  =        temperature  function  defined  by 

equation  (III- 92)  °F 

x  =        axial  coordinate  ft 

y  =        transverse  coordinate  ft 

Y   (x)  -         Bessel  function  of  second  kind  of 

n  j 

order  n 


GREEK  LETTERS 

2  . 
0(.  =        fluid  thermal  diffusivity  ft    /hr 

2 

o(  =        pipe  wall  thermal  diffusivity  -     ft    /hr 

£  =  /V  (used  in  Section  IV)  hr/ft 

f[  =        x/v  (used  in  Section  IV)  hr 

Q  =        fluid  temperature  F 

o 

Q  *  =         radially  averaged  fluid  temperature  F 

J  =  (1-t)/fe      (used  in  Section  IV)  ft 

f  =  fluid  density 

X  -  t-t  (used  in  Section  IV)  hr 

d>  =  fluid  temperature  function  at  interface        F 


lbm/ft 


3C  =        function  defined  by  equation  (V_43) 

•ur  =        pipe  wall  temperature  function  at 

interface 


Also,   additional  notations  are  occasionally  introduced  for 
temporary  use  in  the  analytical  developments  which  follow.     In  such 
cases,    the  new  notations  will  be  adequately  defined  in  the  text.     This 
also  implies  that  the  same  symbol  may  have  different  significance  in 
different  portions  of  the  text.     There  should  be  no  confusion  arising 
from  this  practice.     In  particular,    different  notation  is  frequently 
introduced  for  dummy  variables  of  integration.     Also,    the  symbols 

/5   ,         ^    ,  \    ,         ju  ,    etc.  ,   used  as  eigen  values,    do  not  appear 

in  the  foregoing  list. 

We  also  remark  that  equations  are  numbered  serially  in  each 
section.     In  referring  to  equations,    the  serial  number  alone  will 
indicate  the  correspondingly  numbered  equation  in  the  same  section 
as  the  reference;  otherwise  the  section  number  will  also  appear.     Thus, 
in  Section  V,    for  example,    the  reference  (12)  is  to  equation  (V_12) 
and  the  reference  (II-7)  is  to  equation  (II-7),    that  is,    to  equation  7  of 
Section  II.     A  reference  such  as  IID2  is  to  the  text  of  subsection  IID2 
and  not  to  an  equation. 


ABSTRACT 

Transient  temperature  distribution  in  axisymmetric  systems  of 
incompressible  fluid  flow  through  semi -infinite,    externally  insulated 
cylindrical  pipes  is  considered  under  the  assumption  of  constant  uniform 
velocity  and  constant  material  properties.     The  general  problem,    includ- 
ing all  significant  heat  transfer  mechanisms,    except  radiation,    is 
reduced  to  an  integral  equation  containing  only  one  unknown  function, 
namely  an  interface  temperature  as  a  function  of  time  and  axial  position. 
By  neglecting  one  or  more  heat  transfer  mechanisms,    several  other 
formulations  are  obtained.     These  are  analyzed  so  as  to  provide 
approximate  and  closed  form  solutions  not  previously  available. 
Similar  formulations  and  results  are  obtained  for  the  case  of  fluid 
between  parallel  plates. 
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I.     INTRODUCTION 

A.      PRELIMINARY  REMARKS 

The  problem  of  transient  temperature  fields  in  systems  of  fluid 
flow  through  circular  pipes  or  through  channels  formed  by  parallel 
plates  is  of  considerable  practical  and  theoretical  importance.     On  the 
practical  side,    one  may  consider,    for  example,    the  problem  of  the 
calculation  of  thermal  stresses  in  main  coolant  piping  systems    of 
nuclear  power  plants  in  which  abrupt  changes  in  power  output  may 
cause  very  rapid  variation  in  the  temperature  of  the  coolant  fluid 
leaving  the  reactor.     In  such  cases  one  may  be  concerned  with  the  time 
history  of  the  coolant  temperature  at  some  critical  point,    such  as  a 
pump  or  valve,    downstream  from  the  reactor. 

From  a  theoretical  point  of  view,    the  problems  implied  by  such 
pipe  flows  are  of  great  interest.     They  involve  two  dependent  variables, 
the  temperatures  of  the  fluid  and  of  the  pipe  wall,    each  of  which  is  a 
function  of  three  independent  variables:     time  and  two  space  coordinates. 
These  functions  must  satisfy  appropriate  partial  differential  equations 
and  external  boundary  and  initial  conditions  and  must  also  conform  to 
each  other  in  an  appropriate  way  at  their  common  interface.     Such 
problems  are  called  "conjugated.  "     Few  solutions  for  such  problems 
are  given  in  the  literature.     Those  which  were  previously  known  usually 
involve  drastic  simplifications  such  as  considering  only  steady  state 
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behavior  or  restricting  space  variation  to  one  dimension,    or  considering 
the  pipe  to  have  infinite  thickness,    or  neglecting  important  mechanisms 
of  heat  transfer,    or  a  combination  of  these. 

Accordingly  it  is  important  to  consider  such  problems  in  as  gen- 
eralized a  form  as  possible  both  for  the  purpose  of  obtaining  compact 
alternate  formulations  and  also  for  the  purpose  of  obtaining  closed 
form  and  approximate  solutions  to  problems  in  which  the  degree  of 
simplification  is  less  than  that  involved  in  previously  known  solutions. 
It  is  believed  that  the  work  reported  herein  is  successful  in  both 
regards. 

B.   STATEMENT  OF  THE  GENERAL  PROBLEM 

It  is  appropriate  at  this  point  to  make  an  explicit  statement  of  the 
general  problem  which  will  be  considered  herein.     We  are  concerned 
with  the  flow  of  an  incompressible  fluid  axially  down  the  length  of  a 
semi-infinite  cylindrical  pipe  whose  cross  section  is  a  circular  annulus 
of  constant  finite  wall  thickness.     It  is  presumed  that  axisymmetrical 
conditions  prevail  so  that  the  space  variables  are  radius    r   and  axial 
position  x   from  the  accessible  end  of  the  pipe.     The  outside  of  the 
pipe  is  presumed  to  be  perfectly  insulated  so  that  there  is  no  transfer 
of  heat  through  this  surface.     (This  condition  is  approximately  satisfied 
in  many  cases  of  technical  importance.  ) 

It  is  assumed  that  the  physical  properties  of  both  the  fluid  and  the 
pipe  wall  are  constant  --  that  is,    they  do  not  vary  with  respect  to  time, 
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temperature,    or  position.     We  presume  that  the  thermal  conductivity 
of  the  materials  involved  is  orthotropic;  that  is,   we  provide  for 
different  values  of  conductivity  in  the  radial  and  axial  directions.     This 
permits  later  specialization  in  different  applications,    either  to  the  case 
of  isotropic  conductivity  or  to  the  cases  where  one  or  more  conducting 
mechanisms  are  suppressed.     Subject  to  this  introduction  of  orthotropic 
conductivity,    it  is  assumed  that  Fourier's  law  of  heat  conduction  holds 
in  each  medium.     In  most  of  the  work  done  herein,    we  provide  for  an 
interface  surface  heat  transfer  coefficient  to  account  for  surface 
resistance  or  other  effects.     This  coefficient  may  be  taken  as  infinite 
if  there  is  no  such  resistance.     However,    in  all  cases  it  is  assumed  to 
be  a  constant  which  does  not  depend  on  time,    position,    or  temperature. 

We  have  assumed  the  nature  of  the  fluid  flow,    namely,   we  have 
assumed  that  the  velocity  of  each  element  of  fluid  is  the  same,    a  con- 
stant velocity  V  in  the  axial  direction  x.     This  is  frequently  called 
slug  flow.     It  is  recognized  that  this  is  a  particular  case  which  is  at 
variance  with  the  flow  conditions  which  would  actually  obtain  under 
any  circumstances.     However,    it  is  not  far  removed  from  the  flow 
conditions  which  would  prevail  with  fully  turbulent  flow  and  it  does 
provide  a  necessary  simplification  which  has  permitted  us  to  obtain 
useful  results. 

We  have  confined  attention  to  one  case  of  temperature  variation 
which  can  be  described  as  follows.  It  is  presumed  that  the  fluid  has 
been  entering  the  pipe  at  zero  (or  reference)  temperature  for  a 
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solution  for  the  case  of  cylindrical  geometry.     In  fact,    in  some  cases, 
when  certain  simplifying  assumptions  and  procedures  are  introduced, 
the  two  cases  may  become  one  and  the  same. 

C.      SURVEY  OF  LITERATURE 

The  study  of  heat  transfer  problems  in  systems  of  fluid  flow 
through  pipes  dates  back  to  1883  when  one  of  the  most  celebrated 
problems  of  heat  transfer  was  solved  by  L.    Graetz  [11].  *    In  Graetz's 
problem  the  velocity  field  is  taken  to  be  known,    and  heat  transfer 
takes  place  under  steady  state  conditions  with  pipe  wall  temperature 
distribution  specified  a  priori.     In  recent  years  other  problems  of 
fluid  flow  have  been  treated  by  many  authors.     These  studies  include 
such  effects  as  heat  transfer  by  axial  conduction  and  mechanical  heat 
generation  due  to  friction,    i.  e.  ,  viscous  dissipation.     There  is  an 
abundance  of  literature  in  this  area  of  heat  transfer  engineering. 
However,    these  problems  are  concerned  solely  with  behavior  of  the 
fluid;  the  action  of  the  containment  is  considered  only  in  the  specifica- 
tion of  boundary  conditions  for  the  fluid. 

However,    the  transient  case  of  temperature  field  in  systems  of 
flow  through  pipes  or  ducts  when  the  wall  thermal  state  is  not  specified 
a  priori,   and  it  is  required  that  its  effect  on  the  temperature  field  in 
the  fluid  be  accounted  for,    clearly  fallsin  the  class  of  conjugated 


♦Numbers  in  square  brackets  refer  to  the  Bibliography,    page   154 
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problems  of  heat  transfer.     These  conjugated  problems  are  more 
involved  inasmuch  as  the  mutual  thermal  effects  of  the  wall  and  the 
fluid  must  be  accounted  for.     Accordingly  their  treatment  presents 
serious  mathematical  difficulties.     Even  when  the  velocity  field  is 
assumed  to  be  known,    the  thermal  problems  are  still  of  great  difficulty, 
mainly  because  of  the  coupling  that  is  introduced  by  the  condition  at 
the  common  boundary  of  the  fluid  and  the  wall. 

At  present  there  are  numerous  publications  dealing  with  conjugated 
problems  of  heat  transfer.     In  all  cases,    however,    they  seem  to  be 
particularized  and  as  yet  no  general   method    has  been  obtained  for 
treatment  of  these  problems. 

An  early  treatment  of  conjugate  problems  of  heat  transfer  was 
given  by  T.    L.    Perelman  [25]  who  considered  steady  state  heat  transfer 
to  a  fluid  flowing  around  a  body  containing  internal  heat  sources.     He 
also  considered  the  asymptotic  solution  to  the  type  of  integral  equation 
that  occurs  in  the  analysis  of  such  conjugated  problems:     namely,    a 
Fredholm  equation  of       second  kind  with  singular  kernel  and  infinite 
range.     He  considered  slip  flow  around  a  body  and  also  considered 
the  laminar  boundary  layer  flow  over  a  thin  plate.     Perelman  neglected 
the  axial  heat  conduction  in  the  fluid. 

N.   A.   Avdonin  [4]  considered  transient  behavior  in  a  similar  case 
but  with  no  heat  sources  present.     He  also  neglected  the  axial  conduc- 
tion in  the  fluid.     He  solved  the  problem  by  using  Green's  functions 
and  Laplace  transforms.     He  treated  the  conjugated  problem  by 
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reducing  it  to  a  Volterra  integral  equation  which  he  solved  using 
Laplace  transforms. 

D.    E.    James  and  N.    G.   William  [17]  considered  the  steady  state 
problems  of  Poiseuille  -  Couette  flow  between  parallel  planes.     In 
this  work  the  boundary  planes  were  taken  to  be  of  finite  length  in  the 
direction  of  flow  and  of  finite  thickness  perpendicular  to  the  direction 
of  flow.     These  authors  assumed  an  arbitrary  temperature  distribution 
at  the  fluid -wall  interface  and  solved  the  energy  equations  separately 
for  each  material  for  a  constant  interface  temperature  distribution. 
Then  they  used  Duhamel's  superposition  integral  to  obtain  the  general 
solutions  in  terms  of  the  arbitrary  interface  temperature  distribution. 
Next  applying  the  continuity  condition  for  the  heat  flux  they  obtained  an 
integro-differential  equation  for  the  interface  temperature  distribution. 
They  solved  this  by  the  method  of  collocation. 

Recently  A.    V.    Luikov,    V.    A.   Aleksashenko,    and  A.    A. 
Aleksashenko  [22],    considered  conjugated  steady  state  heat  transfer 
in  circular  and  plane  tubes  for  the  case  of  developed  Poiseuille  velocity 
distribution  allowing  for  mechanical  energy  dissipation.     They  state 
that  general  methods  have  been  developed  for  both  steady  and  unsteady 
heat  transfer,   but  they  give  no  example  or  discussion  for  other  than 
steady  state  cases.     They  essentially  reduce  the  problem  to  a  singular 
integral  equation  for  the  unknown  temperature  at  the  fluid-wall 
interface.     This  is  a  very  interesting  paper. 
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E.    M.   Sparrow  and  F.   N.    DeFarias  [33]  investigated  the  laminar 
flow  through  parallel  plate  channels  under  the  assumption  of  sinu- 
soidally  varying  fluid  inlet  temperature.     This  work  is  in  the  spirit  of 
the  new  results  to  be  reported  herein  inasmuch  as  a  solution  for 
sinus oidally  varying  end  temperature  can  be  related  to  that  for  step 
change  in  end  temperature  and  conversely.     Accordingly  it  is  important 
to  delineate  clearly  what  simplifying  assumptions  they  employed.      They 
neglected  axial  conduction  in  both  fluid  and  pipe  and  they  assume  pipe 
temperature  varies  with  axial  position  but  not  with  transverse  position. 
They  considered  constant  fluid  velocity  as  we  do  herein.     It  happens 
that  these  simplifying  assumptions  do  not  coincide  with  any  set 
treated  in  this  thesis. 

V.   S.   Arpaci  [3]  considered  a  transient  problem  of  the  type  in 
which  we  are  interested.     Assuming  as  we  do,    the  case  of  slug  flow, 
he  uses  the  term  "rod"  to  describe  what  we  call  the  fluid.     He  shifts 
viewpoint  in  an  obvious  and  nonessential  way  by  considering  the  pipe 
to  be  moving  at  constant  velocity  with  its  contents  at  zero  velocity. 
Axial  conductivity  is  neglected  in  both  substances.     Radial  conductivity 
is   accounted  for  in  the  fluid.         Radial  conductivity  of  the  pipe  wall 
is  taken  as  infinite.     The  transient  effect  is  due  to  sudden  and  uniform 
heat  generation  in  the  contents.      Like  ours,    his  pipe  has  perfect 
external  insulation.     He  used  the  method  of  finite  Hankel  transforms 
to  reduce  the  problem  to  the  solution  of  an  integro-differential  equation 
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which  was  solved  by  Laplace  transforms.      For  large  values  of  time  an 
asymptotic  behavior  of  the  temperature  function  was  found,    again  using 
Laplace  transforms. 

Other  significant  contributions  to  the  study  of  conjugated  heat 
transfer  problems  have  been  made  by  M.    D.    Kelleher  and  K.    T.    Yang 
[19],    R.    Siegel  [34],   and  R.   Siegel  and  J.    M.   Savino  [35]. 

F.    E.    Moreland  [23]  employed  numerical  inversion  of  Laplace 
transformation  to  deal  with  heat  exchanger  problems  which  are  similar 
to  our  problems,    except  for  suppression   of  all  heat  diffusion  in  the  fluid. 
D.   D.    Hiep  [13]  extended  this  work  significantly  by  including  the  effect 
of  axial  heat  diffusion.      (His  treatment  is  for  porous  media,   but  our 
formulation  could  be  similarly  generalized  by  use  of  some  fictitious 
material  and  geometrical  properties.  )    His  study  contains  a  large 
bibliography.     The  problem  considered  in  Section  VI  hereof  is  a  special 
case  of  that  considered  by  Hiep.     We  are  able  to  obtain  analytical 
results  whereas  Hiep's  analysis  is  purely  numerical. 

Two  earlier  solutions  which  provide  significant  points  of  departure 
should  also  be  mentioned,    even  though  idealizing  assumptions  made 
the  analyses  very  simple.     One  of  these  is  by  W.    Munk  [24]  who 
assumes  identical  constant  temperature  in  both  pipe  and  contents  at 
any  cross  section.     He  neglects  all  diffusion  of  heat  and  arrives  at  a 
solution  which  repeats  the  time-temperature  history  at  the  entrance 
point,   but  with  a  propagation  velocity  less  than  the  velocity  of  the 
fluid.     The  second,    to  which  many  investigators  have  made  contributions, 
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is  discussed  in  section  15.  3  of  [7]  where  many  references  are  cited. 
Arpaci  [2],    page  353  et  seq.  ,    gives  a  very  compact  treatment  of  this 
problem  in  which  there  is  no  heat  diffusion  in  either  the  axial  or  radial 
direction. 

D.      SUMMARY  AND  OUTLINE  OF  SUCCEEDING  SECTIONS 

The  following  telegraphically  brief  statements  are  intended  to 
provide  a  summary  of  the  content  of  the  succeeding  sections. 

Section  II  -  Derives  basic  equations  which  will  be  employed  in 
subsequent  sections. 

Section  III  -  Considers  radial  and  axial  heat  diffusion  in  each 
material.     Arrives  at  integral  equation  (111-70)  for  slab  geometry  and 
integral  equation  (III- 181)  for  cylindrical  geometry.     No  solution  of 
these  equations  is  attempted. 

Section  IV  _  Radially  averaged  temperatures  are  considered  in 
each  material.     An  ordinary  perturbation  solution  is  constructed 
which  is  valid  "behind  the  (temperature)  front.  "    A  singular  perturba- 
tion solution  is  developed  which  gives  variation  ahead  of  the  front.     A 
scheme  for  matching  these  at  the  front  is  given. 

Section  V  -  By  averaging  the  fluid  temperature  in  the  transverse 

direction  and  neglecting  axial  conduction  in  the  pipe  wall^closed  form 

solutions  are  obtained  in  the  case  of  slab  geometry:     equations  (V-39, 

40).      The  case  of  cylindrical  geometry  turns  out  to  differ  only  as 

regards  the  definition  of  a  certain  auxiliary  function  so  that  the 
solution  has  the  same  expression  as  for  the  case  of  slab  geometry. 
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Section  VI  -  By  taking  radial  averages  for  both  fluid  and  solid 
temperatures  and  by  neglecting  axial  conduction  in  the  solid,    a 
simplified  version  of  the  preceding  problem  is  obtained.      (This  may 
also  be  obtained  by  neglecting  axial  conduction  in  the  solid  in  the  case 
given  in  section  IV).     A  closed  form  solution  is  given:     equation  (VI-21) 
and  (VI-23). 

Section  VII  _  Here  the  problem  is  greatly  simplified  by  assuming 
(in  effect)  that  there  is  no  radial  variation  in  temperature.     Axial 
conduction  in  both  materials  is  considered,    however;     (i.  e.  ,    we  take 
average  fluid  temperature  equal  to  average  solid  temperature  at  any 
section  and  employ  the  equations  of  section  IV).     The  resulting 
solution,    (VII- 10)  is  not  new  but  is  included  here  for  completeness  and 
comparison  purposes. 

Section  VIII  -  This  section  includes  discussion  and  recommendation 
for  further  studv. 
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II.  FUNDAMENTAL  EQUATIONS 

A.   GENERAL  REMARKS 

All  the  analysis  presented  in  this  thesis  relates  to  a  definite 
physical  problem  which  will  be  described  shortly.      The  real  physical 
problem  is  so  complicated  as  to  preclude  the  possibility  of  dealing 
with  it  without  making  certain  idealizations.     In  what  follows  immediately, 
the  problem  will  be  described  in  mathematical  terms  and  the  idealiza- 
tions which  are  common  to  all  the  subsequent  analysis  herein  will  be 
delineated.     The  simplified  problem     which  results  is  still  of  great 
complexity.     It  will  be  further  discussed  in  section  III.     No  solution  of 
this  problem  will  be  presented  but  it  will  be  shown  that  the  problem 
can  be  cast   in  a  more  compact  manner  which  gives  some  promise  of 
being  amenable  to  analysis. 

In  order  to  obtain  additional  results,    it  has  been  found  necessary 
to  make  additional  idealizations  or  simplifications.     These  lead  to 
several  simplified  problems  which  are  specifically  delineated  and 
dealt  with  in  sections  IV  through  VII.      For  all  of  this  work,    however, 
the  derivation  given  immediately  below  and  the  equations  which  result 
therefrom  are  fundamental. 

1.       Physical  Model 

The  basic  physical  problem  of  interest  is  that  of  heat  transfer 
taking  place  between  a  conduit  and  the  fluid  flowing  therein.      The    former 
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is  perfectly  insulated  on  its  exterior  surface.     It  may  be  either  a 
circular  pipe  with  uniform  wall  thickness  or  a  duct  composed  of  two 
parallel  plates.      (In  what  follows  these  two  cases  will  be  referred  to 
as  "cylindrical  geometry"  and  "slab  geometry,  "  respectively.     It  is 
assumed  that  the  fluid  entering  the  conduit  has  had  constant  temperature 
for  sufficient  time  that  the  conduit  and  its  fluid  contents  have  reached 
this  constant  temperature  throughout.     Then  the  temperature  of  the 
fluid  entering  the  conduit  is   suddenly  changed  by  a  finite  amount  and 
is  thereafter  maintained  at  the  new  temperature.     (The  end  surfaces 
of  the  conduit  also  experience  the  same  step  change  in  temperature.  ) 
It  is  required  to  predict  the  unsteady  temperature  fields  in  both  the 
fluid  and  the  conduit  for  all  locations  at  any  subsequent  time. 

Figure  1  illustrates  the  two  geometries  which  are  considered. 
That  for  the  cylindrical  geometry  is  self-explanatory,    except  to  remark 
that  there  is  no  variation  whatsoever  with  respect  to  angular  coordinate. 
The  case  of  slab  geometry  involves  identical  parallel  plates  which  are 
s  emi- infinite  in  extent  in  the  x  direction,    which  is  the  direction  of 
fluid  flow, and  are  infinite  in  the  z  direction.     It  is  assumed  that  there 
is  no  variation  whatsoever  with  respect  to  the  coordinate  z     of  any 
quantity.     In  making  a  mathematical  formulation  of  the  problem,    it 
will  be  convenient  to  take  advantage  of  the  symmetry  with  respect  to 
the  plane  y  =  0;  for  this  reason  the  auxiliary  coordinates  y     and  y 
are  used  as  indicated. 
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To   X  =  oo 


a  -  cylindrical    geometry 


b  -  slab  geometry 


"b  X=<x> 


Figure    I:     Geometry  and  Dimensions  of  Conduits 
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2.       Assumptions 

a.  The  fluid  is  considered  to  be  incompressible. 

b.  The  fluid  is  assumed  to  be  flowing  with  constant  velocity  V 
(slug  flow). 

c.  Material  properties  are  taken  to  be  constant  and  independ- 
ent of  time,    position,    and  temperature. 

d.  It  is  assumed  that  there  may  be  an  additional  thermal 
resistance  (possibly  due  to  a  very  thin  layer  of  foreign  material  at  the 
interface)  the  heat  transfer  coefficient  of  which  is  a  constant. 

e.  Viscous  dissipation  is  assumed  negligible. 

f.  Heat  transfer  by  radiation  is  considered  to  be  negligible. 

g.  Fourier's  law  of  heat  conduction  is  taken  to  be  valid 
throughout. 

h.       The  diffusivity  appearing  in  Fourier's  law  for  fluid  is 
assumed  to  include  not  only  the  molecular  diffusivity  but  also  an  eddy 
diffusivity  which  is  effective.      The  combined  diffusivity  is  assumed 
to  satisfy  assumption  c. 

i.        There  is  variation  only  with  respect  to  time,    the  axial 
dimension,    and  a  radial  dimension,     Figure  la^or  a  transverse 
direction,  Figure   lb. 
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B.      DIFFERENTIAL  EQUATION 
1.       Cylindrical  Geometry 

There  is  no  difficulty  in  establishing,    by  customary  methods 
of  derivation,    the  following  partial  differential  equations:     (1)  which 
describes  the  fluid  and  (2)  which  describes  the  solid  pipe  wall.     We 
have  anticipated  simplifications  which  will  be  made  in  later  sections 
hereof  by  considering  artificial  media  in  which  thermal  conductivities 
in  the  axial  and  transverse  directions  may  differ.      These  conductivities 
are  distinguished  by  a  second  subscript,    a  to  denote  axial  or  t  to  denote 
transverse.     If  we  are  considering  isotropic  materials  (as  we  do  in 
section  III  hereof,    for  example)  we  take  the  axial  conductivity  to  be 
equal  to  the  transverse  conductivity.     However,    if  in  order  to  make 
mathematical  headway  in  a  certain  analysis  we  neglect  axial  conductivity 
in  the  solid  (say),    we  can  obtain  the  corresponding  mathematical  des- 
cription simply  by  setting  K  =0. 

ma 

There  should  be  no  difficulty  in  recognizing  the  significance 
of  the  various  terms  which  appear  in  equations  (1)  and  (2)  below. 

*falfr  it  V  *x  +Kf*~drVTr)  "HCfM  (1) 


k    il+K     ±A(r^I)  =  PC    51  (2) 
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One  may  observe  that  (2)  can  be  obtained  from  (1)  by  taking 
V  =  0,    replacing  9  by  T,    and  replacing  subscript  f  by  subscript  m. 
2.       Slab  Geometry 

The  partial  differential  equations  governing  the  temperature 
distribution  in  the  flowing  fluid  and  in  the  solid  duct  wall  are 


;,    £jg  +|<    ££  -PC  V^  =  f>c   °? 
fa  frx2-     *  ^J|l    h  *     *x     k  S   dt 


for  the  fluid,   and 


£C,  +  kw .  ^:,  ^p^^1"  (   (4) 


ma 


for  the  pipe  wall.     Observe  the  use  of  subscripts   1  and  2  for  the 
coordinate  y,    consistent  with  figure  lb. 

C.      INITIAL  AND  BOUNDARY  CONDITIONS 
1.       Cylindrical  Geometry 
a.       For  Fluid 

(1)  Initial  Condition,  t  =  0.  The  initial  temperature  of 
the  moving  fluid  is  considered  to  be  uniform,  and  without  any  loss  of 
generality,    it  may  be  considered  as  zero. 

0  =  o  o  *  x     <    co        ,         o^rff-Q-         ,       t  =  o 

(2)  Boundary  Condition  at  x  =  0.     At  this  boundary  the 
entering  fluid  temperature  is  given  a  sudden  step  change  to  a  higher 
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value  thereafter.     Without  loss  of  generality,    the  amount  of  increase 
may  be  taken  to  be  unity. 

0  =  1  x-°      j       o  £r  <a     j       t?o 

(3)  Boundary  Conditions  as  x — *oc>.     As  x  approaches 
infinity  the  fluid  temperature  decreases  to  its  initial  value. 

9  — >0  x— >«>,         o<f<a       j        t>o 

(usually,    this  condition  is  employed  in  the  form:     9  remains  bounded 
as  x  — >oo  •  ) 

(4)  Boundary  Condition  at  r  =  0.     It  is  required  that  as 
r — >  0  the  fluid  temperature  remain  finite.     This  condition  may  also 
be  stated  by  requiring  that  the  temperature  distribution  in  the  fluid 
to  be  symmetric  about  the  center  line. 

Either 

or 

-M.    =o  r  =  o     ,      o  <x  <<»        j         +>0 

b.       For  the  Pipe  Wall 

(1)     Initial  Condition,    t  =  0.      The  initial  temperature  of 
the  pipe  wall  is  considered  to  be  uniform  and  equal  to  the  fluid  tem- 
perature at  the  initial  time. 


7"=  o  o<x  <  oo     j       a<<*  <  t> 


t~o 
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(2)  Boundary  Condition  at  x  =  0.     The  pipe  wall  is 
assumed  to  experience  the  same  unit  step  change  of  temperature  as 
that  of  the  fluid. 

T  =  1         x=o        ,       a^r^b         ,         t^o 

(3)  Boundary  Condition  as  x  — *-oo.     Similarly  it  is 
required  that  the  pipe  wall  temperature  approach  its  initial  temperature 
as  x  approaches  infinity. 

T — »o       x  — »■  oo,         a^r      <  t>      ,        t  >  o 

(4)  Boundary  Condition  at  r  =  b.  Since  it  is  assumed 
that  the  pipe  wall  is  perfectly  insulated  at  the  outside  surface,  then 
the  Fourier  heat  conduction  law  requires  that 


vr=o 


O  <  X    <  oo 


r  =  b        ,    t  >  o 


c.       For  the  Interface  Between  Fluid  and  Wall  (r  =  a) 

Continuity  of  heat  flux  is  required  at  the  interface  between 

fluid  and  pipe  wall. 

Contemplating  a  constant  heat  transfer  coefficient  h  at 

this  boundary,   we  may  thus  write 


K     .     ~b&         =  h  (  0  -  T) 
Kmt   S  =  h  (  9  -  T) 


7>r 


o^x<co,     r=a      ,     t^o 
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It  is  well  to  note  at  this  point  that  if  the  heat  transfer  coefficient  h 
approaches  infinity,    then  the  above  conditions  reduce  to  equality  of 
temperature  and  heat  flux,    in  the  following  manner. 


e  =  t 


o  <x  <oo     ,       r*  a 


t  > 


£° 


■ft    dr  **t   ^r 

2.       Slab  Geometry 

The  details  and  physical  significance  of  the  choice  of  initial 
and  boundary  conditions  for  this  case  are  similar  to  those  discussed 
for  the  case  of  cylindrical  geometry.     We  will  therefore  state  them 
mathematically  as  follows. 

a.       For  the  Fluid 

(1)  Initial  Condition  t  =  0. 

9-0  0  <X   <*>  ,        o  <  u    <  <x        j         t>0 

(2)  Boundary  Condition  at  x  =  0. 

6=1  X  =o  j         »<#  <(*       J         *£  o 

(3)  Boundary  Condition  as  x — +00  . 

9  =  °  x— **oo      ?       0<u  <a     j       t>o 

(4)  Boundary  Condition  at  y,   =  a. 


*3. 


s  0 


oix  <oo 


3, 


*  a 


t>0 
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(Cf.   figure  lb;  the  symmetry  about  the  face  y     =  a  provides  this 
condition.  ) 

b.       For  the  Wall 

(1)  Initial  Condition,    t  =  0. 

T  =  0  o  <  x   <  oo         ,       o<y2    |b,       t  =  o 

(2)  Boundary  Condition  at  x  =  0. 

T  =  l  x  =  0  ,       0  -|  y2  «g  b       ,        t  >  o 

(3)  Boundary  Condition  as  x ►co. 

T  =  0  x  — ►  co         ,         °^Y2=b'       t>  O 

(4)  Boundary  Condition  at  y?   =  b. 


=  o 


O    <  X     <  CO 


Y2  =  b      ,       t>  o 


c.       For  the  Interface  Between  Fluid  and  Wall  y.    =  y2  =  0 


■K      .    AL    =  h  (6  -   T) 
Kft      dj>     =  h  (0  -  T) 

*a. 


o  ^  x  <  oo      ,     y  -,  =  Y2  -  o  ,   t  ?  o 


Observe  that  the  algebraic  signs  given  here  are  a  con- 
sequence of  the  fact  that  y.   and  y?  are  opposite  in  direction;  Cf. 
Figure  lb.     Again  it  is  easy  to  observe  that  as  h  approaches  infinity 
the  interface  condition  above  may  be  written  as 

0      =  T 


>      o  ^  x  <co      ,     yj  =  y2  =  o    ,    t  ^  o 


mt  *3i       ft  *3l 
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D.      SPECIALIZATION  OF  ENERGY  EQUATIONS 

In  what  follows  particularizations  of  equations  (1)  and  (2)  for  the 
cylindrical  geometry  and  equations  (3)  and  (4)  for  the  slab  geometry 
are  introduced,    which  involve  various  techniques  for  simplifying  the 
energy  equations.     These  methods  involve  radially  averaging  the  fluid 
and  the  conduit  wall  temperatures  by  integrating  with  respect  to  radius 
(for  cylindrical  geometry)  and  with  respect  to  the  transverse  co- 
ordinates (for  slab  geometry).     Thus  appropriate  definitions  of  average 
temperatures  and  surface  film  heat  transfer  coefficient  are  introduced. 
Further  simplifying  assumptions  are  made,    such  as  neglecting  the 
axial  conduction  in  the  conduit  wall;  these  simplifications  give  rise  to 
interesting  problems.     Also  since  the  difference  between  cylindrical 
and  slab  geometry  disappears  when  both  the  fluid  and  the  wall  tem- 
peratures are  averaged,    only  the  cylindrical  geometry  is  considered 
in  this  case. 

1.       Cylindrical  Geometry 

a.       Radial  Averaging  of  Temperatures 

We  consider  isotropic  conductivity  in  both  media: 


K     "   K„*=:    Km  (6) 

We  also  define  the  average  temperatures 
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a 


r 


and 

I 


I  redr 


9    =    -?- --     -±-    I  redr  (7) 


r 


(8) 


£ 


and  we  introduce  an  "equivalent"  surface  film  coefficient 


(9) 


Note  that  h*  relates  heat  flux  to  difference  between  average 
fluid  temperature  and  average  wall  temperature. 

Then,    by  integrating  equations   (1)  and  (2)  with  respect  to 
radius,   between  appropriate  limits,    and  introducing  the  definitions 
above,   we  find 

*»-^M   -    «f«  *Yt  (10) 

for  the  fluid  mean  temperature,   and 

fir"  „_■_  2I*_  £f  (e»_r*)  ,ii, 
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for  the  pipe  wall  mean  temperature.     Equations  (10)  and  (11)  are 
essentially  those  considered  by  Hiep  [13];  although  his  equations 
appear  to  provide  for  more  generality  (in  order  to  apply  to  porous 
media),   a  change  of  variable  shows  that  they  are  mathematically  the 
same.     The  introduction  of  h     may  legitimately  be  objected  to  on  the 
grounds  that  physically  it  would  be  quite  difficult  to  evaluate.      However, 
its  use  provides  a  degree  of  generalization  compared  to  the  physically 
reasonable  simplifying  assumption  that  radial  variation  of  temperature 
in  the  fluid  is  so  small  that  9''"  =  0  (r  =  a)  =  T  (r  =  a).     Moreover, 
even  in  cases  where  h — *co,    the  use  of  h'  <.  co    permits  taking  into 
account  some  resistance  to  radial  heat  flow  in  the  fluid. 

The  initial  and  boundary  conditions  for  (10)  and  (11)  follow 
naturally  from  those  of  (1)  and  (2). 

(1)  Initial  Condition  t  =  0. 

6"  =  T*  =  o  o  ^  x  <oo  ,  t  =  o 

(2)  Boundary  Condition  at  x  =  0. 

9*  =  T*  =1  x  =  o  ,  t  >  o 

(3)  Boundary  Condition  as  x  — >oo. 

e':— >0     ,     T* — >o      x  — »-co  ,  tt° 

b.       Case  where  h'~  Approaches  Infinity 

An  interesting  special  case  may  be  obtained  if  it  is  assumed 
that  radial  diffusion  is  so  great  that  the  temperature  is  essentially 
constant  throughout  any  particular  cross  section.     This  is  accomplished 
by  letting  h"  approach  infinity.      Then  from  (11)  we  obtain 
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(e*_T*)  =  *=£=  (±?L-£r  )  (12) 


Substituting  for  (0     -  T":)  from  (12)  above  in  (10),    and  noting  that  as  h''~ 
approaches  infinity,    then  6  =  T    ,    yields 

for  the  common  fluid  and  conduit  wall  temperature  distribution.  The 
initial  and  boundary  conditions  for  (13)  follow  naturally  from  those  of 
(10)  and  (11). 

(1)  Initial  Condition,    t  =  0. 

T*  =  0  o  <  X   <oo      ,  t  =  0 

(2)  Boundary  Condition  at  x  =  0. 

T*  =  1  x  =  0      ,         t  >  0 

(3)  Boundary  Condition  as  x  — >ao. 

T*  — >   0  x  — >oo  ,         t  >  0 

The  solution  of  this  system  is  treated  in  section  VII. 
c.       Case  where  K       =  0 

By  neglecting  axial  conduction  in  the  solid  pipe  wall,    a 
tractable  system  is  obtained.     We  rewrite  equation  (10)  with  no  change 

tie  -*  "*£*«  JL  U\  •*!?  (6>*-T*)  do) 
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Multiplying  equation  (11)  by    Km   and  allowing      K     to  approach  zero, 
yields 


.«     -* 


Tit         P  ^    c 


(14) 


The  initial  and  boundary  conditions  are 

(1)  Initial  Condition  at  t  =  0. 

0*  =  T*  =  0  O^X<oo      j       t=.o 

(2)  Boundary  Condition  at  x  -  0. 

9*  =  1  X=o         ^        t£  o 

(3)  Boundary  Condition  as  x  — *oo  • 

9  — v  O  y — +00       j       t^° 

The  solution  of  this  system  is  dealt  with  in  section  VI  . 
d.      Another  Special  Case 

If  radial  averaging  is  done  only  for  the  fluid  and  variation 
with  respect  to  radius  is  retained  in  the  solid  pipe  wall  an  interesting 
problem  results.     However,    this  situation  leads  to  an  integro- 
differential  equation  the  importance  of  which  does  not  seem  to  be 
commensurate  with  its  difficulty.     Accordingly  an  additional  simplifica. 
tion  is  made:    axial  conduction  is  neglected  in  the  solid  pipe  wall.     To 
accomplish  this  we  set    K         =  0    in  equation  (2)  and  we  also  write 
Kmt  =  Km.      For  the  fluid  we  write  Kfa  =  Kft  =  Kf.     Note  that  the 
radial  effect  is  absorbed  in  the  radial  averaging  for  the  fluid.     It  is 
necessary  to  define  a  new  equivalent  heat  transfer  coefficient  fr  as 
follows. 
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-Km2l   ,-K*6    -!(••- T)         rf     r-a 


}r  f    * 


Then  equations  (1)  and  (2)  with  0'    still  defined  by  (7)  reduce  to 


(is; 


>*-(tU) 


(16) 


and 


il(4),il  (17, 

(17)  is  exactly  the  form  taken  by  (2)  with  Kma  =0.      (16)  is  the  same 
as   10)  except  for  the  new  definition  of  surface  film  coefficient.     The 
equivalent  heat  transfer  coefficient     h     ,    defined  by  (15)  relates  inter- 
face heat  flux  to  the  difference  between  conduit  surface  temperature 
and  fluid  mean  temperature;  compare  with  h'~  defined  by  (9)-      Clearly 
h   V  h    .     The  reasons  for  introducing    h    are  covered  in  the  remarks 
following  (11). 

The  appropriate  boundary  and  initial  conditions  are: 

(1)     Initial  Condition.      The  initial  conditions  for  equations 
(16)  and  (17)  are  exactly  the  same  as  those  for  equations  (1)  and  (2) 

0*  =  0  O  <X   <  00        ;        t*  ° 


T  =  0 
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(2)  Boundary  Condition  at  x  =  0.      The  same  conditions  as 
stated  for  equations  (1)  still  apply  to  equation  (16): 

0*  =  1  X=0  J        *  >° 

We  observe  that  since  equation  (17)  is  explicitly 

independent  of  x,    it  is  therefore  not  required  to  specify  any  condition 

at  x  =  0  for  equation  (17). 

(3)  Boundary  Condition  as  x_xjo.     It  is  again  required 
that  Q^  remain  finite;  viz. 

0*  =  0  X  — +0O     )         t  >o 

(4)  Boundary  Condition  at  r  =  b.     Assuming  perfect 
insulation  at  the  outer  surface  of  the  cylinder,    then 

"fcr 

(5)  Boundary  Condition  at  r  =  a.     At  the  interface  of 
fluid  and  pipe  wall  solid  the  continuity  of  heat  flux  requires 

or 

ft 

We  note  here  that  if    h    approaches  infinity,    then 

the  condition  of  continuity  of  heat  flux  above  reduces  to  equality  of 
temperature  at  the  fluid  and  wall  interface. 

This  system  is  analyzed  in  VB1.     See  discussion  at 
end  of  IID2  for  explanation  of  why  the  system  of  IID2  is  treated  first. 
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2.       Slab  Geometry 

For  all  the  simplified  cases  in  Section  IID1  in  which  average 
temperatures  in  both  fluid  and  pipe  wall  are  considered,    the  geometry 
of  the  situation  is  fully  described  by  the  parameters  A      A     ,    and  P. 
Thus,    specializing  the  case  of  slab  geometry  in  the  same  ways  leads, 
as  one  may  easily  verify,    to  exactly  the  same  equations.     However, 
the  case  treated  in  IIDld  still  retains  radial  variation  in  the  pipe  wall. 
Accordingly,   we  expect  that  the  slab  geometry  counterpart  will  be 
described  by  slightly  different  equations,   which  we  now  proceed  to 
obtain. 

a.       Differential  Equation 

Treating  equations  (3)  and  (4)  for  the  case  of  flat  wall 


duct  similarly  we  define 

a 

a*-  ^       i  fed, 


0 
Then  considering  isotropic  conductivity  in  the  fluid, 


integration  of  (3)  yields 


111     vj*£  -  JL  It  +  jL  (V  t) 


(19) 


for  the  fluid.     Similarly  with  the  assumption  of  zero  axial  conductivity 
in  the  wall,    we  obtain 
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ti  =  j.  21  ,20) 

for  the  solid  wall.     Note  the  introduction  of  h  in  (19);  see  remarks 
following  (17). 

b.  Fluid;  Initial  and  Boundary  Conditions 

(1)  Initial  Condition  t  =  0. 

9*  =  0  0  ^  x    <  co  ,         t  =  0 

(2)  Boundary  Conditions  at  x  =  0. 

9=1  x  =  0  ,t^0 

(3)  Boundary  Conditions  as  x — ►  co  . 

e"~ >  0  x  — >.  co  ,         t  >  0 

c.  Solid  Wall;  Initial  and  Boundary  Conditions. 

(1)  Initial  Condition,    t  =  0. 

T  =  0  0  5;  y2  <f.  b  ,         t  =  0 

(2)  Boundary  Condition  at       2  =  b. 

^T  =  0  Y0=b  ,         t  >  0 

*&.  2 

(3)  Boundary  Condition  at  ^2  =  0. 

■a- 

-K         "ftT     =    h    (9    "    T)    y   =  0         ,         t>  0 
m  ^ 

dy2 

This  system  which  is  similar  to  that  obtained  in  subsection  IIDld  is 
somewhat  easier  to  deal  with.     Accordingly,    it  is  treated  first,    in  VA 
while  that  for  the  system  IIDld  is  dealt  with  in  VB. 


40 


III.     INTEGRAL  EQUATION  FORMULATION 

A.      INTRODUCTION  AND  PLAN  OF  ANALYSIS 

Equations  (II-3)  and  (II-4),    with  Kfa  =  Kft  =  Kf  and  Kma  =  Kmt  =  Km 
so  as  to  provide  isotropic  thermal  conductivity,    involve  two  dependent 
variables,    0(x,  y,  t)  and  T(x,  y,t),    each  involving  three  independent 
variables,    along  with  a  formidable  number  of  boundary  and  initial 
conditions.     Similar  remarks  hold  for  the  case  of  cylindrical  geometry. 
These  are  complicated  systems  to  deal  with  and  we  have  achieved 
solutions  which  are  useful  for  calculating  purposes  only  by  making 
various  simplifying  assumptions.      These  solutions  are  given  in  other 
sections  of  this  thesis.     In  this  section  we  will  represent  the  system 
of  equations   (II-3)  and  (II-4)  in  another  form  which  involves  only  one 
unknown  function  of  two  independent  variables.     No  additional  conditions 
are  required  for  the  complete  specification  of  this  unknown  function. 
The  form  of  this  equation,    (70),    is  that  of  an  integral  equation  of  what 
is  known  as   "first  type"  (that  is,    the  unknown  function  does  not  appear 
outside  an  integral  sign)  with  complicated  kernels  and  forcing  functions. 

Although  the  method  of  procedure  is  essentially  the  same  for  the 
case  of  cylindrical  geometry  as  for  the  case  of  slab  geometry,    the 
details  are  sufficiently  different  to  justify  treating  the  two  cases 
separately.     We  treat  the  case  of  slab  geometry  first  and  in  subsection 
HID  we  return  to  the  case  of  cylindrical  geometry.     In  each  case  the 
method  of  procedure  is  essentially  as  described  in  the  next  paragraph. 
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We  define  an  unknown  function    a)  (x,  t)  which  represents  the 
(unknown)  temperature  of  the  fluid  at  the  interface  and  a  similar 
function  "UV  (x,  t)  which  represents  the  unknown  temperature  of  the 
solid  at  the  interface.      For  the  fluid  we  solve  the  problem  of  finding 
the  interior  temperature  9(x,  y,  t)  when  the  surface  temperature  is 
given  by   yCx,  t),    subject  to  other  appropriate  boundary  and  initial 
conditions.     It  is  convenient  to  consider  this  solution  as  the  sum  of 
two  solutions,    one  involving  a  homogeneous  boundary  condition  at  the 
interface  and  accommodating  the  boundary  condition  at  the  end,   x  =  0, 
and  a  second  solution  involving  a  homogeneous  boundary  condition 
for  x  =  0  and  the  unknown  function  ^(x,  t)  at  the  interface.     The 
solution  to  the  latter  problem  is  expressed  in  terms  of  an  integral 
involving  the  function <b(x,  t);  it  is  obtained  by  use  of  a  procedure 
which  is  due  to  A.    N.    Lowan  [20,21].     We  proceed  in  the  same  way 
to  express  the  solution  for   the  solid  temperature  T(x,  y,  t)  in  terms 
of  the  unknown  surface  temperature   "^(x,  t).      The  interface  heat 
balance  condition  provides  two  relations,    one  of  which  may  be  used 
to  express  the  function  ll/  (x,  t)  in  terms  of  the  function      <f>  (x,  t). 
The  other  relation  then  becomes  the  integral  equation  for  the  unknown 
function    <p  (x,  t). 

B.      AUXILIARY  PROBLEMS  (SLAB  GEOMETRY) 

We  will  deal  with  the  problem  of  the  temperature  distribution 
in  the  fluid  first.      Then,    by  taking  the  flow  velocity  as  zero  and  by 
using  subscript  m  (for  the  solid)  rather  than  f  (for  the  fluid),     and 
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by  changing  the  designation  of  the  surface  temperature  distribution 
from    <fr  (x,  t)  to      y(x,t),   we  can  very  quickly  obtain  the  corresponding 
solution  for  the  temperature  distribution  in  the  solid.     We  immediately 
write 

6(x,  y,t)  =  uj(x,  y,t)  +  u2(x,y,  t)  (1) 

where  u-(x,  y,t)  are  the  solutions  to  the  two  auxiliary  problems 
considered  in  the  remainder  of  this  subsection. 

1.       Fluid;  Homogeneous  Boundary  Condition  at  Interface 

For  simplicity  in  this  subsection  we  will  omit  subscripts, 
and  will  replace  them  when  we  combine  solutions:  o(  will  denote 

0^£,     u    will  denote  uj,    y  will  denote  y,    (Cf.    Figure   lb).      Thus  the 
problem  of  interest  is 


ty*      bx2-  ~x    2>x     '  «    It 


(2) 


The  boundary  and  initial  conditions  are 


u   =    1 

x  =  0 

(Ky^a 

t>  0 

(3) 

u  — 0 

X  ->oo 

Oly^a 

t£  0 

(4) 

u  =  0 

O^x  <  oo 

y  =  0 

t>  0 

(5) 

*9 

0*x  <  oo 

y  =  a 

t>0 

(6) 

u  =  0 

CKx  <oo 

0<y#  a 

t  =  0 

(7) 
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Using  an  overbar  to  denote  Laplace  transform  with  respect  to  time, 
we  obtain 


M  +!5    _J^^t  =_i-a  (8) 


This  incorporates  the  initial  condition.     All  the  boundary  conditions 

will  be  the  same  in  Laplace  transform  form  except  that  at    x  =  0 

which  becomes  (jL  =  _-     Thus  assuming  a  solution  of  the  form 

S 

—  Ax 

U  -  A  e      '       sin(yiy)  (9) 

we  observe  that  all  conditions  are  satisfied  except  that  for    x  =  0, 
and  the  condition  of  zero  slope  at    y  =  a.     Imposing  the  latter  yields 
the  characteristic  values 

U     =  (n+l/2)lT/a  n  =  0,1,2,...  (10) 

Because  of  linearity  we  can  write 


An  e      *U(K%)  (11) 

n*o 


Substitution  of  (11)  into  (8)  yields 


x„=.v.  i^ft^| 


2<< 


(12) 


The  negative  root  is  discarded  in  order  to  assure  boundedness  as 

x  — ».  oo  •      The  condition  at  x  =  0  and  use  of  the  orthogonality  property 

of  trignometric  functions  yields 
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A     =      -I 
n  S 


Substituting  from  (13)  and  (12)  into  (11)  we  obtain 


(13) 


oo 


»lrO 


(14) 


where 


*1  1  v        2 


(15) 


The  inversion  of  (14)  may  be  obtained  using  standard  tables 

[29].     Thus 

Co 

U(x,y,t)  =       >        Y,(y)  X^x.t) 
where 


(x,y,  t)  =      V     Y^y)  X^x, 


n-o 


(17) 


x,(x,o« 


1      2o< 

—  e 


e^(^#jiU(^J 


(18) 


In  these  expressions,    although  there  is  explicit  dependence  upon  the 
integer  n,    this  is  not  reflected  in  the  notation  on  the  left.      This  is 
done  merely  for  simplicity  of  notation  and  will  also  be  done  with 
similar  definitions  to  follow.     In  every  case  the  symbols  are  used 
under  a  sign  of  summation  with  respect  to  n. 
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For  future  purposes  we  evaluate  the  partial  derivative  of 
(16)  with  respect  to  y  and  evaluate  at  y  =  0. 


oo 


o        n=o 


where 


y'(0)  -  1  (20) 

As  the  reader  may  come  to  appreciate,   when  the  many  preliminary 
results  are  assembled  into  the  final  equation,    a  compact  notation  is 
convenient.     Functions  Y,  and  X,  and  constants      u      are  the  first  of 
several  such  notations  which  will  be  introduced. 

2.       Fluid;  Non-homogeneous   Boundary  Condition  at  Interface 
In  what  follows    U  stands  for  U?(x,  y,  t)  and  it  is  governed  by 
the  same  equation  (2).     However,    the  following  changes  are  made  in 
the  boundary  conditions  (3)  and  (5) 

U=0  (rather  than  u  =   1)        x  =  0        0<y|a      t>  0 

U=  <f>(x,  t)    (rather  than  u  =  0)     O^x^oo        y=0  t^O 

where      $>  (x,  t)         =M  <oo      for  all  x  and  t.     We  further  make  the 
following  substitution  in  order  to  drop  out  the  convection  term. 

u(x,y,t)  -  W(x,y,t)  exp  f ¥*  .  V±  }  (21) 

Then  equation  (2)  and  the  appropriate  boundary  and  initial  conditions 
reduce  to  the  following  equation  and  conditions  for  W(x,y,t). 
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dx*         bj2 


*  at 


(22) 


W  =  0 


Xso        ^     °£c)^a       ;         *£° 


(23) 


W<jo 


X-^co    j     °<%i<>-     ,         t>o  (24) 


Wrg(x,t)  ^X<«>        ,  J^o  ,  t|<) 


(25) 


|^=o  o<'X<oo     ,  y=a    ,?  t>0  (26) 


W=  o 


o<)C<a>     ^      o<3  <a    ^  t=, 


(27) 


where 


g(x,tj  =^(x,t;  ^xp 


r  z, 

- vx  +  it 


(28) 


Again  taking  Laplace  transformations  with  respect  to  t  and  accounting 
for  the  initial  condition  we  obtain 


bx2-       Stf        * 


(29) 


All  the  boundary  conditions  will  be  exactly  the  same  as  above  except 
(25),    which  becomes     W    =  g  (x,  s). 
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Following  the  method  of  solution  used  by  Lowan  [20,  21]  we 
consider  the  function 

oo  oo 

W^S)  *JffJn*jV*)j}%JS)'£*l*X  A-^d^  (30) 

oo 


where 


F(a,u;jJg   *tff[A(«-B)3  (31) 


and 

X2=yu2+-|.  (32) 

Note  that  M  and  t   are  dummy  variables  of  integration.     Direct 

substitution  shows  that  (29)  is  satisfied.     It  is  easy  to  verify  that  all 

the  boundary  and  initial  conditions  are  satisfied.     In  particular  we 

note  that 

CO  oo 

o  ° 

by  the  Fourier  identity.     As   regards  boundedness 

CO   OO 

'  o  o  " 

This  completes  the  proof  that  (30)  is  the  desired  solution  of  (29) 
which  satisfies  all  the  auxiliary  conditions. 
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Continuing  to  follow  the  procedure  described  by  Lowan  [20,21], 
we  define 


CO 

-St 


/**  £(**>*;* s  f<a'^x)  -  rbvffil)    (33: 


i.  e.  ,    F„     (a,  y;t)  is  the  inverse  of  F   (a,  y;  A  ).     The  subscript     U 
indicates  the  dependence  upon  the  parameter    u  ;  the  dependence  of  F 
upon  s  is  made  clear  by  expressing    A  in  terms  of  s,   as  above.     The 
convolution  theorem  then  gives 

oo  oo  t 

O    O  o 

Further  defining 


(35) 


then 


W*>*>*>  .  H*.l,t)  ,36, 


Thus  substituting  from  (36)  in  (34)  yields 
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oo  CO  t 


If  we  now  take  the  Laplace  transform  of  (35),   we  obtain 


*     iCs) 

where 


Y(5J  =   cos-k  [X(a-gj]  ;         E(sj=  ca5*(AaJ 


(39) 


The  singularities  of  (38)  are  simple  poles  corresponding  to 
the  zeros  of  cosh(   X  a);  °f  course,    these  are  functions  of  the  variable 
u    .     Note  that  cosh  (A  a)  =   1  +     A  a/2!  +    \  a4/4  !  +  .  .  .   which  converge; 

everywhere  in  the  finite  plane.      Thus,    although   ^\     itself  involves  a 

2 

square  root      ,         \     does  not,   and  there  is  no  difficulty  to  be  en- 
countered dealing  with  branches.     Thus,    the  residue  theorem  affords 
the  desired  inverse  of  (3  8) 


*<•'»*»•£♦£ 


oo  Snt 

«*^£-  (40) 


s  r s? 


5=S, 
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The  summation  in  (40)  extends  over  the  roots  of 

cosh  (A  a)  =  0  (41) 

namely  the  numbers 


x»-^l?  *  (n+»/2)^r=<>, 


«-*,l,2,...  (42) 


(Do  not  confuse  the  constants      Mn.  of  (10)  with  the  continuous  variable 
/>  of  (30)).     Thus 


(f'+O 


(43) 


where  Sn  are  the  poles  of  (38)  in  excess  of  the  simple  pole  at  S  =  0. 
Note  that  they  are  functions  of  the  variable     V    .     Therefore  the  formal 
inverse  of (38)  is 

The  following  derivative  is  needed  in  (37) 

-  H    **■  =  ^  >  (-0  m,  c<*  |//nC^; I  e  (45) 
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Hence,    substituting  from  (45)  in  (37)  we  obtain 


oo 


H« 


W(x^t)  =  ^  Vc-OV    cos  [M  a -a  )]fT 

L — i  a    o 

«  =  o 

-      OO         2 


00    t         i 
o  ~o 


(46) 


(-"PI    . 


•  > 


Jnd 


1d^ 


In  order  to  complete  the  analysis,    substitute  (28)  in  (46)  and  then 
substitute  the  result  into  (21).      Thus  we  obtain 

t — i  o  0 


t(5#H) 


(47) 


.3  J 


We  are  able  to  carry  out  the  integration  with  respect  to   U      so  that  we 
finally  obtain 


Q3,  OO    t 

n=o         °  ° 


where 


(48) 


^ 


e/*. 


1  wn 


i} 


-  e^c  I 


(49) 


(50) 
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For  later  purposes  we  need 

V 


f   ~bU 


which  is   evaluated  below. 


>"  V- 


du 

* 


OO        Co   x 
J  M  =  o    '       ■*     a    a 


(51) 


In  order  to  assure  that  there  is  no  misinterpretation,    we 
remark  that  the  integration  with  respect  to  ^    is  from  zero  to  infinity 
and  that  with  respect  to    v^    extends  from  zero  to  t. 
3.       Combination  of  Solutions  for  the  Fluid 

Simple  addition  of  (16)  and  (48)  yields 


«»  oo  t 


(52) 


Likewise  by  superimposing  (19)  and  (51)  we  obtain, 


00  eo  t 


(53) 


where,    having  restored  the  appropriate  subscripts, 


Y,0,)=  ^^ 


X,(x,t)  = 


Vx 


+ 


« 


XU)  =  ^  ^ 


Ma-a) 


(54) 


(55) 


(56) 
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.J 


& 


exp 


-  exp 


HI 


(57) 


Also,    of  course,    (10)  and  (15)  define    n    and      u 

"n  rn 

It  is  well  to  pause  at  this  point  to  review  the  physical  nature 
of  the  problem  for  which  (52)  provides  the  solution.      This  is  given 
below.      The  differential  equation  is 


(H-3) 


The  auxiliary  conditions  are 

X->oo     ^     °  =  3,^a 

D=    O  0<    X   <00      ,        0  <  #,   <<2 

Solution  for  the  Solid  Wall 


0- 

i 

e< 

00 

©- 

4>(x 

,t> 

0 

4. 


,       t>o 


The  solution  for  the  similar  problem  for  the  solid  wall  is 
easily  obtained  from  the  preceding  solution  for  the  fluid.     All  that  is 
necessary  to  do  is  to  write  T  in  place  of  9,  y  in  place  of  4*    ,    replace 
subscript  f  by  subscript         m,    replace  a  by  b,    take  V  =  0,    and 
introduce  a  new  symbol     ^n  (defined  below)  in  place  of    MM  . 
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T(*JM*t' 


(58: 


where 


Y^(^)  = 


(59) 


xz(*,o=  4- 


(60) 


Y(,]a  (-0Weos[xn(fc-»j] 


(61) 


X2     /3"    """^"l      I 


-(fcx) 


-<xp 


:62 


X«-    (n+l/i)£ 


(63) 


(Do  not  confuse       Xn  above  with       A  n  used  temporarily  in  (42)) 


«*       _  I 


3T 

*& 


oi"        — '  o  o 

»*0 


(64) 
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C.      DERIVATION  OF  INTEGRAL  EQUATION  (SLAB  GEOMETRY) 

Equations  (52)  and  (58)  thereby  give  the  fluid  and  pipe  wall 
temperature  distributions  respectively.     However  these  solutions  are 
in  terms  of  the  arbitrary  interface  temperature  functions       T  (x,  t) 
and        "Y  (x,  t).     We  now  return  to  the  actual  boundary  conditions  that 
prevail  at  the  interface,    viz. 


Kffa^(t-Y) 


8  =  °      >     °$*  <°°      )    t^°     (65) 


and 


%\r°       }     otx<oo     >     t>o    (66; 


(Note:     The  signs  on  the  left  hand  side  of  (65)  and  (66)  are  consistent 
with  the  choice  of  coordinates  used  in  Figure  lb).     Solving  for  Y  from 
(65)  in  terms  of     <t>    ,   we  obtain 


»•♦-*«&■ 


(67) 


Substitution  from  (53)  in  (67)  yields 


'  oo  oo  t 


(68) 


Ki(x*1>M)d,ld5 
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Furthermore,    substitution  of  (68)  in  (66)  yields 


-K 


vr 


9) 


Substitution  of  y  (x,  t)  from  (68)  into  (64)  and  the  result  into  (69) 
yields 


oo 


£ 


»J=o 


O  < 


CO 


Xi(»jt)+nKx(K^*>n)' 


0     O 


f(l.*l)-$)     X,feM) 


m-=t> 


tjf*(«,n-fi)  &(?,*,  t-i^M^* 


O    0 


H<±\  + 


(7o: 


It  should  be  clear  that  we  could,  with  equal  facility,  express  the 
problem  in  terms  of  "U/  (x,  t)  instead  of  <t>  (x,  t)  as  above.  The 
distinction  is  trivial. 

Equation  (70)  is  the  desired  integral  equation,    the  only  unknown 
in  which  is  the  function    <£>    (x,t)  which  represents  the  surface  tem- 
perature of  the  fluid. 

Equation  (70)  is  considerably  simplified  if  there  is  no  surface 
resistance.     Letting  h,    the  heat  transfer  coefficient,    approach 
infinity,    reduces   (70)  to 
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'  f  o    o 

n  =  o 

(71) 


[K.^Mh  —  ^^'i) 


Jad^ 


for  h  =  co  . 

D.      AUXILIARY  PROBLEMS  (CYLINDRICAL  GEOMETRY) 

For  the  case  of  cylindrical  geometry,    the  treatment  is  quite 
similar  to  that  which  has  preceded.     However,    it  is  not  possible  to 
derive  the  equations  for  the  pipe  wall  by  simple  manipulation  of  the 
equation  for  the  fluid  as  was  done  above.     Thus  it  is  appropriate  to 
consider  four  auxiliary  problems   similar  to  those  dealt  with  before, 
two  for  the  fluid,    arriving  at  equation  (165),   and  two  for  the  solid, 
arriving  at  equation  (170).      These  equations  will  then  finally  be 
combined  so  as  to  arrive  at  the  desired  integral  equation.      By  suitably 
redefining  the  several  functions  which  are  involved,    it  will  be  possible 
to  express  this  equation,    (181)  exactly  in  the  form  of  (70). 

1 .       Fluid;  Homogeneous  Boundary  Condition  at  Interfa c e 
Consider  the  problem  mathematically  stated  as 


lxz     <*    dx       r    3r        ?>rx  ~  *    dt  (/  ' 
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with  the  boundary  and  initial  conditions: 

U  -  1  X=o  j      o<r  ia.     >         *><>  (73) 


U=  o 


X-*«>         *    aS   C  ia      j  t>o  (74) 


—  =  0  o<K<oo  ,  r  so      j  t>o  (75) 


*r 


U.-    o 


«<x  <<*>       5        r^a   }       t>o       (76) 


u-=  o  o<x<oo       _,     *<r<a    }        t*o       (77) 

Assuming  a  solution  of  the  form 

K   -     AJ.(^)   X|(K,tJ  (78) 

and  substituting  (78)  in  (72)  yields 

<£&  „XM)'  y?2X,  =  JL  Mi  (79) 

7>x*         *      ^  *     bt 

Equation  (78)  already  satisfies  condition  (75).      Condition  (76)  yields 

the  characteristic  equation 

J.(PaJ=0  (80) 

the  roots  of  which  are  denoted  by   n   ,    n=l,  2,  .  .  .       .     Because  of  the 

linearity  we  can  write 

00 

u=  VAnJ.(/?nr)X/(x,t)  (8i) 

The  solution  of  (79)  is  obtained  by  taking  Laplace  transform  with 
respect  to  time,    which  yields  a  second  order  homogeneous  ordinary 
differential  equation  the  solution  of  which  is 
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X,(x,t)  =   B  -**? 


£-*/(£)'+ /»M 


(82) 


In  (82)  the  exponential  with  positive  exponent  has  been  discarded  so 
that  condition  (74)  may  be  satisfied.  Substitution  of  (82)  in  (81)  and 
absorbing  B  into  An  ,    we  obtain 


oo 


f«-*\(U+ti+i 


(83: 


U  =    VAnJe(/3„r]^p 

Thus  at  this  point  we  observe  that  the  inverse  of  (83)  satisfies  (72), 
(74),    (75),    (76)  and  (77).     Condition  (73)  provides  a  means  of  evaluating 
An  .      Taking  x  =  0  (73)  gives 


oo 


1  =    yAnJo(^r) 


(84) 


*-l 


multiplying  by  r  J0  (  /?   Y*  )dr  and  using  the  orthogonality  relation  of 


Bessel  functions; 

a 


fj.{P«r)Jo{P*r)r<ir  =  ±  J,  (fima) 


(85) 


if  n=m  and  zero  if  n^  m,    one  easily  finds 


A„  = 


I       2 


S     P'*     J,(W) 


Thus  we  now  have 


nn 


Mt»*) 


e*p 


%-^nfpR] 


(86) 


(87) 
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The  inversion  of  (87)  is  tabulated  [29];  viz. 


oo 


(88) 


where 


e  "  er^c  ( 2=2^* )  +  e  ""erf  c  f ^^) 


(89) 


90; 


For  future  purposes  we  need  the  partial  derivative  of  (88) 
with  respect  to  radius  evaluated  at  r  =  a. 


6]  =-£x,(x,tj 

'r=  o.        L — • 


(91) 


2.       Fluid;  Non-homogeneous  B.C.  at  Interface 

The  problem  considered  here  is  the  same  as  that  in  the 

preceding  subsection  except  that  the  following  changes  are  made  to 

the  boundary  conditions  (73)  and  (76). 

U  =  0  (rather  than  u  =   1)  x=*o    J     0<r<a     f     t>o 

U  =  4>(x,t)       (rather  than  u  =  0)       0<x<Cb,  r=  a     J     t>c 

where    \<p  (x,  t)  I  <M  <(  00         for  all  x,  t.     A  substitution  of  the  form 


"(r.V)--  W(.,x,OM£-^j 


(92) 
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reduces  (72)  to 

ax2-     r  df     ^r2-    oc   &t 

and  all  the  boundary  and  initial  conditions  will  remain  the  same 
except  that  at  r  =  a  we  have 

W(a,  x,  t)  =  cj>  (X,  t)  exp   (  "^  +  ^j  )        =  g(x,  t)  (94) 

Taking  Laplace  transform  of  (9  3),    we  obtain 

<l$/        J_  W  +  *W    =  S  ^  (95) 

bxl       r  br       dr2-      * 

Again  all  the  boundary  conditions  will  be  the  same  except  (94)  which 
in  Laplace  transforms  will  become 

Wia^XjS)  a  3(*jS)  c?4T^^o>  r=a  (96) 

Following  the  method  of  solution  used  by  Lowan  we  consider 
the  function 

vfajs)* tfffFfavi)  |(?^s;/^/*xAi.^d/*ci|      (97) 

0  o 
where 


(98) 


and 


X  =  ^  +.  x.  (99: 
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Note  that    f*  and    *|  are  dummy  variables  of  integration.     Direct 
substitution  shows  that  (95)  is  satisfied  as  well  as  all  the  boundary- 
conditions.     In  particular  once  again  we  note  that 
oo  oo 

oo  (100) 

=  I<*,s; 

by  the  Fourier  identity.     As   regards  boundedness, 


00    OO 


W(f,X,s)|^    ^fT/W^x  /W^djuJ^   =  l«t 


(101) 


This  completes  the  proof  that  (97)  is  a  solution  of  (95)  which  satisfies 
all  auxiliary  conditions.      Let  us  now  define 


co 


je  t  E  {**;  *)*«  F  (^  r)  \)  =  {far;  {fiT) 


(102) 


i.  e.  ,    F n  (a,  r;t)  is  the  inverse  of  F  (a,  r;     X   ).     Note  again  that  the 
subscript  I*    indicates  the  dependence  upon  the  parameter     P     ;  the 
dependence  of  F    upon  s  is  clearly  indicated  by  expressing    X     in 
terms  of  s,    as  shown  in  (102).     The  convolution  theorem  then  gives 


oo  oo 


Wfax,*)*  ifffopxAUfti  ^J|J^(ff,C7)g(^t-»i)di      (103: 


o   o 
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Further,    defining 


Gpfarst)  mffpi^rstfJ* 


(104) 


gives 


(105) 


Substituting  from  (105)  in  (103)  yields 


0  0  0 


(106 


If  we  now  take  the  Laplace  transform  of  (104),    we  obtain 


oo 


where 


o 

-  i  mi 

s      Us) 


(107) 


Y(0=   L(>r) 


2(0  =     LUaJ  (108) 
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Using  the  (temporary)  substitution    \  =  L  ft     ,    (98)  may  be  written 


F     = 


Joiill 


(109) 


However  J0(§  )  is  an  entire  function  having  only  real  zeros.     Thus 
the  only  singularities  of  F  are  poles  at  the  real  roots,    ft  n  of 

Jo(/W  =  0  (110) 

The  corresponding  values  of  s  are: 

s.,=  -*(A/£)  (111) 

Following  Lowan,    we  write  the  desired  inverse  as 


QA 


We|       "l    bS     JS=S,, 


.  Urj    2 


I0(raJ  "a-L  [pn\ 


pMwe 


r\-\ 


(112) 


Thus  taking  the  derivative  of  (112)  and  substituting  into  (106)  yields. 


-*A 


£>         .  A^n^x^M"?  d/i 


dnd 


(113) 


1d? 


It  is  possible  to  carry  out  the  integration  with  respect  to  p. 
in  (113)  and  obtain 
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rr 


(114) 


where 


R2(r)    = 


JAP*') 


(115) 


KiN>*A)-^-s^  e 


"ST    /3n      I 


.   e 


(ii6; 


i 


e<p 


4*>| 


4*1 


In 


Note  that  K    (x,t ,  t,  *[  )  given  in  (116)  differs  from  K^(x,t ,  t,  \  )  given 
by  (50)  only  in  the  eigenvalues;  otherwise  they  are  identical. 

For  later  purposes  we  need  the  partial  derivative  of  (114) 
with  respect  to  r,    evaluated  at  r  ■  a, 


°o       oo    t 


=- y/y>(?.t-i)K,(x,TE,t,.i)«i*<K    (in) 


3.       Solid  Pipe;  Homogeneous  Boundary  Condition  at  Interface 
Consider  the  problem  stated  mathematically  as 


(118) 


with  the  boundary  and  initial  conditions: 
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U=   i  x=o       ,     a<r  <h         }    -t  >0  (119j 


U  =  o 


x~>co   ^,  a<r<{>  ,    t>o  (120) 


°^X<ooj        r=b  _>t>o  (i2i) 


^=0  0<X<<*       j        f=  a-  3     t**  (122) 

a=©  o<x<oo      j,     a^r   <  i>         ,    ts-o  (123) 

We  first  define  the  following  notational  device  for  a  linear  combination 
of  Bessel  functions  of  the  first  and  second  kind  of  an  integer  order. 
Let 

&{*>*}*)  -JK{^)Yj^)-^r)Jrt(tb)  (124) 

Now  assume  a  solution  of  the  form 

u  =    A£i(r,l>jr)  Xz  (*jt)  (125) 

Substitution  of  (125)  in  (118)  yields 

The  function  (125)  already  satisfies  (121),    i.  e.  , 


67 


"br 


=  *£(M;*)  -o 


(127) 


r=b 


Condition  (122)  yields  the  characteristic  equation 


x6t  [ol,  b;tf)  =0 


(128) 


the  roots  of  which  are  denoted  by      Y„    ,    n  =  1,  2,    .  .  .    .     Because  of 
linearity  we  can  therefore  write 


00 


U     =     V  *«&{?;*>  ITS)  Xl(x,t) 


(129) 


n=\ 


To  solve  (126)  we  take  Laplace  transform  with  respect  to 
time.     This  yields  a  second  order  ordinary  homogeneous  differential 
equation,    the  solution  of  which  is 


X2(>os)=     Bexp  f-x  /tfTs/rf 


(130) 


In  (130)  the  exponential  with  positive  exponent  has  been  discarded  so 

that  condition  (120)  may  be  satisfied.     Substitution  of  (130)  in  (129) 

and  absorbing  B  into  A    ,    we  obtain 
6  n' 


oo 


5(r,X>sJ  =  YAn  fifa *>&)**?     -x/tf+V* 


(131) 


»i=J 


Thus  at  this  point  we  observe  that  the  inverse  of  (131)  satisfies  (118), 
(120),    (121),    (122)  and  (123).      Condition  (119)  provides  a  means  of 
evaluating  A    .     Taking  x  =  0  (119)  gives 
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oo 


\  -    Y  A.  fi  (r.bjy.) 


(132) 


n=i 


It  is  not  too  difficult  to  show  that  (124)  has  the  usual  orthogonality- 
properties  of  Bessel  functions  with  the  weighting  function  r.     Thus 
we  merely  present  the  result  as 

if  n  =  m  and  zero  if  n  ^  m.     Use  of  this  orthogonality  property  yields 


A    =  -L 

M"      s 


.2  ^l 


.OL^faWj-bfifaM-X 


(134) 


Substitution  of  (134)  in  (131)  yields 


u 


=1 


4-^(-x^lJ 


(135) 


Note  that  the  parameter  s  appears  only  in  the  last  term.     The  inverse 


of  (135)  is  tabulated.     Hence, 

oo 


u  (r,x,  t)  =\  R3(r)  X^x,  t) 


where 


»»=t 


(136: 


(137) 


xxcx.*j  =  ^"W  (*^±)  ♦  * V  (*$££)    d3 


Y*F 
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The  partial  derivative  of  (136)  with  respect  to  r  and  evaluated 
at  r  =  a  is 


00 


(139) 


where 


»,(«)  = 


V 


# 


This  result  will  also  be  used  later. 


-1 


(140) 


4.       Solid  Wall;  Non-homogeneous  Boundary  Condition  at  Interface 
We  consider  the  immediately  preceding  problem  with  the 
following  changes  in  boundary  conditions 

U  =  0  (rather  than  u  =   1)      x  =  0         a<  f  <  b  t>  o 

U  =  y  (x,  t)     (rather  than  u  =  0)      x  =  0  r  =  <X  t>o 


w 


her  i 


Y  (  X  ,  t )  I     ^M  ^oo    for  all    x  and  t.     After  Laplace  transforma 


tion  we  have 


.2- 


dx*       r   zr\     lr  '        <* 


(141) 


with  the  conditions 
U  =o 


X  -o 


X-+oo 


a  *r  <  b 


ol  <r  <  t 


(142) 
(143) 
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~    =    o  0<K    <«  ,  r=  t  (144) 

U  =  y(x,s)  o*x<c0       j        r^a  (145) 


The  solution  of  (141)  may  be  constructed  following  the  procedure  used 

by  Lowan.     Thus,    consider  the  function 

00  oo 


o  o 


where 


and 


FfrAr;*;-       I'MWti+«-MU")  (147) 

I*(H   k,(Ab)  +  Ktf(Aaj  J((AbJ 


,2.  i 

A    «  ^    +  S/<*  (148) 


(I     and  K     are  modified  Bessel  functions  of  first  and  second  kind  of 
n  n 

order  n).     Observe  that    u    and    *C     are  dummy  variables  of  integration. 

Direct  substitution  shows  that  (141)  is  satisfied.     It  is  also 
easy  to  verify  that  all  the  boundary  conditions  are  satisfied.     In 
particular  we  note  that 

CO  OO 

U  (a,X,0=  fffff  (?,*)  A^^X  AkjOg  <JH|  (149) 

o  o 
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by  the  Fourier  identity.     As  regards  boundedness, 


CD  co 


0  o 


CO 


(150) 


This  completes  the  proof  that  (146)  is  the  desired  solution  of  (141)  and 
satisfies  all  the  auxiliary  conditions.     Again  we  follow  the  procedure 
indicated  by  Lowan,    and  the  sequence  of  analysis  is  precisely  that 
employed  in  subsection  IIIB2  for  the  case  of  slab  geometry;  the  present 
case  is  different  only  in  details.     Accordingly,   we  abbreviate  the 
treatment  and  write  immediately, 


co  oa  ^ 

o  o  o  *■ 


iH 


151) 


and 


<vM,r;t)  ./"  [l  f(^b,rjXj]=/ 


J.  Yili 
s     2(f) 


152) 


where 


Z(0   -L(xajK,Uy  *   K6(xa)r,(M>} 


(153) 
(154) 


Using  the  temporary  substitution       X=  iY  »    (147)  may  be  written 

Ffe^-y).    JWTtW-JifrbjXW^^f,^        (155) 
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where  we  have  used  the  definition  (124).     It  is  reasonable  to  assume 
that  the  function      J£     (a,b;  If  )  considered  as  a  function  of    t    has 
only  real  zeros.      Thus  the  only  singularities  of  F    are  poles  at  the 
real  roots,      y    ,    of 


>C    (^6jYnJ   »  o 


(156) 


The  corresponding  values  of  s  are, 


S.   -  K  (tf  ...  /j 


(157) 


Following  the  procedure  used  by  Lowan,    for  the  inversion  of  (152), 
we  obtain 


sn-t 


*"'  "  "Si 


s  =  s* 


(158) 


ro(pajK,(pb)  +  K.(HaJI,(pfll) 


c» 


1       *1 


X^MJe*0^* 


(159) 


Thus  taking  the  derivative  of  (15  9)  and  substituting  into  (151)  yields. 


oo 


eo 


"(r.x,t)=  y^WJjYStt-^Kiix^,^)^       (160) 


where 
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RH(r)=         *w*^*  d-f  bj  *""  J 


ks  (*>;*»;+ *£(*,*>;*«) 


(161) 


K»Oo?,Mj«  |f.A 


exp 


■(T-x) 


-«o<.p 


f  (162) 


For  future  use  we  require  the  following 


00      '        ?r 


f-a   n*  i 


(163) 


where 


R«<*)- 


*r. 


>C(ftjt;yJ 


+  a 


5.       Combination  of  Solutions  (Cylindrical  Geometry) 
a.       For  the  Fluid 


(164) 


Simple  superposition  of  the  solutions  (88)  and  (114)  yields 
™  cot 

e(<-,x,t>2^  [R»^)Xi{x,*J  +  RiC«ojf^(f>a)K,(x,t,t,t;dt^     d65) 

where  having  restored  the  appropriate  subscripts 
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Ri(r)=R,(r)  =    _L 


J.(M 


*   J,  (pO 


(166) 


11    f  -^„x 


167) 


(168) 


exp 


Rf-xJ 


2-. 


-  -e*p 


-Jig) 


Similarly   superposing  (91)  and  (117)  we  obtain 


}r   J 


(169) 


b.       For  the  Solid  Pipe 

Again  simple  superposition  of  (136)  and  (160)  yields 
oo  a»  i: 

T(r,x,t)=  ^s(f)X^tJ  +  iWt'^  (HO 

v»_.  o    0 


where  again  restoring  the  appropriate  subscripts 


**    ^>b^.)-bV,(fc,b;^ 


(171) 


n    (f)  a  __  _b_  ^(^bjfrQ 


y"     fc&(a,b;Y.i>aSJ(^M*) 


(172) 


KzM^nJ  =  >hf  "T 


f\ 


fe*-p 


-CF-x/ 


(173) 


-  «*i 


11 
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similarly  superposing  (139)  and  (163)  yields 


oo  »  •» 


where 


R,(aJ  = 


R*(«)» 


1 


Jb      qO|   (bj  b^tfw)      "     -J 

a  .e.t^bjwj" 

V„Z 


(175) 


(176) 


E.      DERIVATION  OF  INTEGRAL  EQUATION  (CYLINDRICAL  GEOMETRY] 

We  now  proceed  as  we  did  in  section   IIIC  for  the  case  of  slab 
geometry.  Equations   (165)  and  (170)  give  the  fluid  and  solid  temperature 
distribution  each  in  terms  of  the  arbitrarily  preassigned  function 
$  (x,  t)  and   "Y"  (x,  t).     If  we  now  return  to  the  actual  boundary  condition, 
namely  that  of  continuity  of  heat  flux  at  the  interface,   we  have 


-K, 


11 


=  -8(<P-yJ 


(177) 


*         • 


r*a 


and 


-K 


>n 


51 


=  *  (f-Y) 


(178) 
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and  solving  (177)  for  y  in  terms  of   <j>  ,    we  obtain. 


■ft 


19 


+  <*> 


(1791 


rod 


Substitution  of  (169)  into  (179)  yields. 

00  cot 


Mai  «    O 


180) 


Further  substitution  for^  from  (180)  above  into  (174)  and  the  result 
into  (178)  yields 


oo    r  <*>   «■ 

L — I  O  O 


(181) 


rn=/ 


oo 


JJ  j,sI  o    © 

The  same  remarks  as  those  following  equation  (70)  for  the  integral 
equation  of  slab  geometry    may  be  made     here  too.     Similarly  therefore 
if  we  let  h  approach  infinity,    we  obtain 


-o 


oo      t  oo  t 


(182; 
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This  is  equivalent  to  (71).     Note  that  (182)  can  be  written  more 
compactly  as 


nf     *  "* 


2__ [x  («,t) -►  Jfifai-ti K*(^j^\)<h^ 


n-i  o  o 


=  o 


(183) 


where 


(184) 


K  (M*M)cK|(>oijM)+r"R/,WK,,(x,5A'0  d85) 


Kf 


A  similar  abbreviation  may  be  made  for  (71). 

The  reader  should  note  that  (70)  and  (181)  are  of  exactly  the 
same  form;  their  only  difference  is  in  the  particular  form  assumed  by 
the  kernel  and  driving  function  and  in  the  fact  that  the  counterparts  of 
RJa)  and  R4(a)  appearing  here  are  simply  unity  in  the  case  of  (70). 

F.      REMARKS  ON  SOLVING  THE  INTEGRAL  EQUATION 

It  is  obvious  that  the  integral  equation  (70)  or  (71)  for  the  case  of 
slab  geometry,    or  (181)  or  (182)  for  the  case  of  cylindrical  geometry 
will  not  be  solved  easily.     It  involves  multiple  integrals,    one  with 
variable  upper  limit  and  one  with  fixed  infinite  upper  limit.      The  kernel 
and  driving  functions,    although  well  behaved,    are  quite  complicated. 
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One  of  the  ways  of  dealing  with  such  a  complicated  problem  is  to  do 
what  has  been  done  in  this  thesis,    namely,    make  additional  assumptions 
so  as  to  simplify  the  problem  to  the  point  that  it  becomes  analytically 
tractable.     Our  attention  has  been  exclusively  directed  in  this  way  and 
we  have  made  no  serious  attempt  to  deal  with  these  fundamental  integral 
equations.     However,   the  following  remarks  and  conjectures  may  be  of 
interest. 

All  the  analysis  herein  shows  that  while    ^>  (x,  t)  is  a  smooth 
function,    almost  all  its  variation  takes  place  in  the  immediate 
neighborhood  of  the  wave  front  x  =  Vt.     It  is  conjectured  that  a  contour 
map  of  <J>  (x»  t)  would  be  similar  to  that  shown  in  Figure  2.     In  this 
figure  the  contours  are  shown  as  straight  lines  merely  because  the 


x  =  Vt 


Figure  2.     Conjectured  Contours  of    c|?  (x,  t). 
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nature  and  degree  of  curvature  are  unknown.      It  is  tempting  to  consider 
replacing  the  boundary  and  initial  conditions    <f>  =   1  along  x  =  0  and 

<f>   =  0  along  t  =  0,    by  the  conditions    <|>  =   1  along  OA  and     <f>  =  0  along 
OB  so  as  to  provide  a  field  in  which  the  "action"  of  <J>  (x,  t)  is  more 
evenly  distributed.      This  might  perhaps  be  useful  in  an  attempt  at 
solution  via  Gallerkin's  (or  a  similar)  method  in  which  a  solution  is 
assumed  which  involves  unknown  scalar  constants  and  these  are 
determined  so  as  to  minimize  an  error  functional.     In  such  a  procedure, 
it  is  evident  that  simplifications  would  result  from  assuming  the  con- 
tour lines  to  be  straight. 
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IV.     METHOD  OF  PERTURBATION 

In  this  section  we  employ  the  method  of  perturbation  to  treat  the 
problem  described  by  equations  (11-10)  and  (11-11)  and  the  associated 
boundary  and  initial  conditions.     The  reader  will  recall  that  these 
equations  resulted  from  having  taken  radial  (or  transverse)  averages 
of  temperatures  in  the  fluid  and  in  the  pipe  wall.     Thus  radial  (or 
transverse)  conductivity  is,    in  effect,    absorbed  into  the  surface  heat 
transfer  coefficient  h*  defined  in  equation  (II- 9).     The  influence  of 
axial  conductivity  appears  explicitly  in  the  following  analysis. 

An  outline  of  what  will  be  accomplished  in  this  section  may  be 
helpful  to  the  reader.      First,    an  ordinary  perturbation  expansion  is 
made  in  terms  of  the  parameter     6  =  1/V.     It  turns  out  that  only  even 
powers  of  €.  are  involved.      The  zeroth  order  solutions  are  identical 
to  well  known  solutions  described  in  section  15-3  of  [7].     The  effect 
of  longitudinal  conductivity  does  not  appear  in  this  zero  order  solution 
but  it  does  appear  in  subsequent  solutions.     It  is  clear  from  the  sequence 
of  higher  order  perturbations  that  the  solutions  fail  at  and  ahead  of 
the  wave  front,    x  =Vt.     Severe  singularities  are  encountered  at  the 
wave  front  and  the  solutions  are  identically  zero  ahead  of  the  wave 
front. 

Thus  a  singular  perturbation  development  is  introduced  to  explain 
the  behavior  ahead  of  the  front.      The  lowest  order  solution  for  this 
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singular  perturbation  expansion  is  similar  to  the  known  solution  for  a 
moving  fluid,    with  diffusion  and  with  an  isothermal  boundary  condition. 
Terms  in  this  singular  perturbation  development  involve  unknown 
functions  of  time  which  are  determined  by  a  process  of  matching  with 
the  solutions  valid  behind  the  front. 

It  should  be  remarked  that  the  system  of  equations  dealt  with  in 
this  section  is  obtainable  from  either  slab  geometry  or  cylindrical 
geometry.     The  process  of  radial  or  transverse  averaging  of  tem- 
peratures erases  the  distinction  between  the  two  geometries. 

A.      ORDINARY  PERTURBATION  ANALYSIS 

Equation  (11-10)  and  (11-11)  together  with  associated  conditions 
are  repeated  here  for  reference  purposes.     The  notation  has  been 
slightly  changed;  we  omit  the  asterisks  which  indicated  radial  averages 
and  we  introduce  the  convenient  constants  b     and  b?  defined  below. 

In  the  fluid,    we  have 

«ffe-v!f  -M^c-t)  (i) 


where 


b,-±£ 


(2) 


subject  to  the  conditions 
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9-0 

For  the  solid  we  have 


X  =    0        i  t  >   0 


X.—+00  j  t^>  o 


<K<oo     j  *t=o 


o( 


&-£-M«-t) 


^x2-        Bt 


where 


b      *P 


PC  A 


m 


subject  to  the  conditions. 


(3) 


(4) 


(5) 


T-i 

TUo 


X-*oo  ,         t  >  o 


o  <x  <oo    ,  -t-  o 


(6) 


It  is  convenient  to  normalize  with  the  following  change  of 
variables 
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t  =  t 

t=  x/y   =  £  x  ; 

so  that  equation  (1)  and  conditions  (3)  become 


e=  */y 


*   "bf-    an    "at    ^     ' 


9  =  1 

9-+o 
9=o 


1=o       ^         t^o 


°£\  1  oo 


t~o 


and  similarly  equation  (4)  and  conditions  (6)  become 


(7) 


(8) 


(9! 


c*  ill  ^^I-bW^-r) 


do; 


r-i 

1=° 

t>a 

r-^o 

i-^a? 

t>o 

T~o 

0  £l  <  oo 

t=o 

(11) 


Equations   (8)  and  (10)  with  conditions  (9)  and  (11)  can  be  approached  by 
the  use  of  Laplace  Transformations  as  was  done  by  Hiep  [13],    but 
there  seems  to  be  no  way  of  solving  the  transformed  equations  or 
performing  the  inversion  other  than  by  strictly  numerical  procedures. 
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This  system  of  equations  may  be  considered  as  depending  on  the 

2 

parameter  £      which  appears  in  (8)  and  (10).      For  large  values  of  V, 

this  parameter  is  small  so  that  it  is  reasonable  to  consider  a  per- 

2 
turbation  expansion  in  £    . 

2 

(Note  that  £  itself  does  not  appear  in  (8)  through  (11),    only   £    ; 

2 
thus  the  expansion  we  develop  is  in  powers  of   £    ,    that  is,    only  even 

powers  of  £  are  involved.      There  are  other  ways  of  showing  that  odd 

powers  should  not  appear  in  this  perturbation  expansion,    but  the  one 

given  here  is  sufficient.  )     Thus, 


H\>tii)s  B*[\»t)+  £  0j(l>tj+...  (12) 


T(^t;^j   ^('btJ^Ti^t)*...  (13) 

2 
Substitution  of  (12)  and  (13)  in  (8)  and  (10)  and  equating  powers  of    £   , 

the  following  sequence  of  equations  are  obtained 


|&+3&+b,  (0.-T,)    =0  (H) 


££«  +  "*£»  +  bi(ai--n)  =(*£§?  (i5 


^r  +  lr  +  b'(ei-"&)  =1  ^rr    -»    £i*    (16) 


for  the  fluid  and 
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25 
it 


-t>L(90-To)    -  o 


2£  -b2  (9,-77)  =  *,^?£ 


for  the  solid.      The  boundary  and  initial  conditions  are 


and 


T.(<w*)  -i 


%=   O       j 

t>o 

<?^1  <©•  j 

fc*P 

V0  j 

0  ^  1  ^°° 

1=«>      ; 

t£o 

0^1<<»      j 

t  =o 

(17) 


(18) 


<->2.         (19) 


(20) 


(21) 


(22) 


["E(i/t)  =  ° 


1=o 


t> 


0  ^1  <  oo   ^  t"  =  o 


6i     (23: 


Since  the  order  of  these  equations  is  less  than  that  of  the  original 
equations,    only  one  boundary  condition  can  be  imposed.     It  is  not 
expected  that  the  solutions  will  be  uniformly  valid.     Starting  with  (14) 
and  (17),    we  take  the  Laplace  transformation  with  respect  to  t.      The 
transformed  equations  are 
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ft  So 

7>n 


+  S0fl  +  b,  (ea-f0)    = 


=  a 


(24) 


(s+biK    =    bz  0< 


(25) 


Expressing  Tain  (25)  in  terms  of    ~§o   ,    substituting  in  (24),    and  solving 
the  first  order  homogeneous  differential  equation  yields 


SStb,  '  LV  S+bi/ 


b.b?. 


(26) 


°.(i>*)  *  j  **? 


bibz. 


-i(^b,-|^-) 


(27) 


The  procedure  for  inversion  of  (26)  and  (27)  is  as  indicated  in 
Example  7-5  of  [2].     The  results  are 


0o  (Vt)   =o 


for       t<^       (28) 


0»  (l,t)  =  e 


e'-I.fc^K^r^tli^ttJ)^ 


for      t^»|       (29) 


and 


1(1.*)  =  ° 


r 


T.U.*).  wiMpb2aio(^)^ 


for      *<*£       (30) 
for        t>1       (31) 
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where  a  =■  b^b^l  and  f  =  t-  v\_  .  Note  the  use  of  a  in  this  sense;  it 
should  not  be  confused  with  a  linear  dimension.  It  may  easily  be 
shown  that       Q0  {*[,  t)  and  To(  *l  ,  t)  both  tend  to  unity  as  t— *■<»  . 

It  is  clear  that  the  two  solutions  obtained  for  T0  (  \  ,  t)  and    Q0  ( *[,  t) 
have  a  wave  like  character  with  no  disturbance  for       \  ^  t.     It  is  well 
to  pause  a  moment  to  refer  back  to  the  physics  of  the  problem.      For 
fixed    n,   (i.  e.  ,    fixed  value  of  x)  no  disturbance  takes  place  until  t    >  v\, 
or  for  fixed  t,    there  is  no  disturbance  for  the  region    >\  >  t.     In  what 
follows  the  condition     X    -   t    -  >\  =  0  is  referred  to  as  the  "wave  front.  " 
The  pipe  wall  temperature  T0  (  *L  ,  t)  is  continuous  at  the  wave  front 
with  discontinuous  slope,    while  the  fluid  temperature     0,  (  *l  ,  t)  is 
itself  discontinuous  at  the  wave  front. 

The  first  corrections      0,  (  *l  t  t)  and  Tj     (  v\  ,  t)  are  governed  by 
equations   similar  to  those  for       Q  (  i\  ,  t)  and  T0  (  Y[  ,  t)  except  that  these 
new  equations  are  nonhomogeneous  with  the  second  derivative  of 

80(  *[  »  t)  and  T0  (  v\  ,  t)  as  forcing  functions.  These  non-homogeneous 
terms  vanish  ahead  of  the  wave  front,    and  thus  no  correction  ahead  of 
the  wave  front  is  expected  in  the  sequence  of  solutions      90  ,    T0  ,     B\, 
Tj,    .  .  .     .     In  the  next  subsection  we  will  demonstrate  this  by  actually 
constructing  the  solutions    Q^  and  T     and  will  see  that  the  present 
perturbation  expansion  fails  at  the  front. 
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B.     HIGHER  ORDER  SOLUTIONS 


The  corrections     Q,  (  ^  ,  t)  and        "T  ( «\  ,  t)  are  also  obtained  by- 
using  Laplace  transformation.      From  (18) 


tf* 


(s+U)T{    =   bid,  +  *w  qjs 


and  from  (26)  we  have 


bVo  _      b; 


ty*       S(S+b2) 


Likewise  from  (15' 


(^-^M^^^iSfe)] 


M  +  (*+b,)5, 


=    bjTi  +  oC 


tfe0 


*  W 


and  from  (27) 


LVk  S  +  bi' 


(32) 


(33) 


(34) 


(35) 


Substituting  (33)  in  (32)  and  solving  for  T,   in  terms  of    Q(    and  then 
substituting  the  result  in  (34)  results  in  a  first  order  non-homogeneous 
differential  equation  for     0,    (  *]  ,  s). 


f^+b'-s»=S 


K^W 


-*l(^b'-^) 


(36) 


The  solution  of  (36)  is  easy  to  obtain. 
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e,=s-n 


c* 


(-^■t)^h(s+b'-^) 


(37) 


It  is  now  merely  a  question  of  substitution  of      0  .   from  (37)  and 
we  have 


from  (33)  in  (32)  to  obtain   T^(   Y^   ,  s).     After  rearranging 


TS  bi 
i  =  ; 


(S  +  biJ 


2  r 

v      Stbi'  lv  sW. 


(38) 


To     invert  (37)  it  is  convenient  to  express    8  ,(    ^  ,s)  as  follows 


9,  =  s  e   .  e    .  £ 


a     4 


where 


"■£*.(* +«».Jn 


n=o 


(39) 


=   2«  1*1*1 


k2        =     °^*lb'  *  ^mblbi^ 

k3     =    2^mbtb2»a 


(40) 


Similarly 


=z  -H,   -n$      s 


«     s 


>j=i 


(41) 
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where 


f   l 
=    2L  *  \i\  bj  \  +  oCm  bi 


k3*       =     b,b2(o<^b,  -H2^m4-  ^mbZ^) 


(42) 


k   * 


=    °i^bU 


For  the  inversion    it         is  useful  to  employ  the  formula 


/ 


-/ 


s  -e 


n-i 


=  (£)2I.(^J 


(43) 


for  integer  n   >   1.      For  n  =  0,    the  formula  is  made  valid  by  adding 
the  term      8  (t)  to  the  right  side;  this  is  the  Dirac  delta  function.     In 
order  to  express  the  result  in  terms  of  modified  Bessel  function  of 
first  kind  order  zero  and  one  only,   we  use  the  following  recursion 
formula 

In+1(x)  =  In_!(x)  -  C2n/x)  In(x)  (44) 

for  modified    Bessel  functions  (Cf.W.    G.    Bickley  [5],    Number  (20a)). 
It  is  also  found  convenient  to  introduce  the  function 


h(  /,n;t)  =  e"bt(a/t)       I   (2  fit) 

n       T 

for  which  the  following  differention  formulas  are  valid. 

dh(l,n;t)    =h(/+1/2,n-l;t)  _(b  +  ^/t+n/2t)  h(  £  ,  n;t) 
dt 

=h(  ^  +1/Z,  n+l;t)  -(b+  fyt-n/2t)  h(  i  ,  n;t) 

In  order  to  invert  (39)  consider  the  following  function 


(45) 


(46) 
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f(n>s)=  e)insn 


-i 


F(v*J-i    fM 


(47) 


*=o 


using  the  formulas, 


i 


-I 


r  d    1 


=  S(tJ  +{£  U(i{&) 


(48) 


and 


/ 


O-n 


-n    f 

S.e 


w-i 


-(£)    l    L,(2^"0  «>=' 


(49) 


yields, 


F(-i.tja„sw+U(|)""rH(2^j        .so, 


yi=»o 


Next  using  the  recursion  formula  (44),    (50)  may  be  expressed  as 


F(*l,t)=kQ    &  (t)  +m0  l0{2y[at   Hir^jFl,  (z{oi)  (51) 


where 


m0     =    kx  +  (Va)  (k 


2k 


[m,     =     \  o^  (a/t)  +  k2  +    2  Ji|_  -    Ji3_  +k<//a) 


(52) 


Thus  we  have 


i 


-/ 


a     ^ 


-bit 

-  e 


&(*;+moI0(z^at )  +  (53) 
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Now  let 


£(%S  +  bz)    »S-  f  (l,S«.b») 


Thus 


(54) 


-*- 


f(*+j 


>-/ 


+  £    sf(l>s  +  b2)-F(o+) 


■*S(t)  +  jtF(i,t) 


To  perform  the  indicated  differentiation,    (46)  is  used;  thus 
-fc»t 


6(i»t;-  c 


(55) 


where 


o"lS(tJ  t^l(2k,-bi)g(tj  +  C„I„(l{ai)  +     (56) 

ClIl(2.{S) 


57) 


ti-  -(V^0(«/*A  ^i(2b--b^(a/0'/z  + 

1/ 


3  ,i 


(«»«M)(t/a) 


V2 


It  is  now  easy  to  obtain  the  inversion  of  (39),    which  may  be  written 


as 
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-1b'    -V     - 


(58) 


Thus 


9iK,t)=  o 

e,(l,t)  =  e    G(n^t-i) 


for  t  <  ^ 

for  t>  1 


(59) 


proceeding  in  a  similar  way,    consider  the  function 

JL  S 


f*(l'S) 


rtsi 


Thus 


Fft>t)=]fc(£]VL,-.(*fS*) 


«=» 


(60) 


(61) 


Expanding  the  sum  in  (61)  and  using  the  recursion  formula  (44)  yields 
F*h»*)*  ™*    I.  (Z^)  +  /I    «ri|(Z^*J  (62) 


where 


Wo 


*      2*4        6i 


,z  „* 


*♦(*)(•«.-£-$•  Ml)*. 


Again  let 


(63) 


5*(v«+fc0  =  sf*u>s+k0 


(64) 
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Thus 


-I 


-* 


6    (l.J+klJ 


r'r  * 


(65) 


=  i  S(t)  +5F*(i»tJ 


The  indicated  differentiation  in  (65)  is  obtained  using  (46).      Thus 


G  h,t)  =  e 


-b2t 


*  ^(O  +C*Ie{z{a±)  +  crii(2|5EJ 


(66) 


where 


v-o 


=  (^tz-*-  2b|t>2^^«o<   b2  1  J  +  (t/aJ  *Mb,bz  • 


i        .3,2, 


ft/a)o<wblb2rl 


(67) 


K  '/ 

C,    =o<b2^(|j2+    oCm(bIb\i_bl/fl-6b1bl/^-2b2/n-Vrl^2bJ;b2J 


+   ^b.b^  (b,-2bz+l/a)(t/<2j 


l/z 


+ 


(t^^l)^!  -2IML-3/,,)  (t/fl) 


?/2 


It  is  again  easy  to  obtain  the  inversion  of  (41),    which  may  be  written  as 


(68) 
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Thus 


Ti(tt.t)  -     e     .6  (i,t-»i) 


for  t"4^ 


for  ^^1 


(69) 


Although  the  inversion  of      Q,    (  ^  ,  s)  and  T,(    1  ,  s)  were  performed  in 
as  compact  a  manner  as  possible,    it  is  desirable  to  rewrite      0|  (  *\    ,  t) 
and  Tj(  v\   ,  t)  as 


ftK^J-TIll'tJ    =° 


for        t  <  1|  (70 


-bi*i  _b2.tr  r         / 


Colo(i(orj  +Cili(2jf5r) 


for        t>r|_  (71) 


and 


-ibi  -b2r 


li(vt)=e  .e 


oc^bzgw  +  c^^/orj+c^/afj 


for        t>1 


(72) 


It  is  of  interest  for  future  discussion  to  obtain  the  limiting  values  of 
9^  n  ,  t)  and  T    (  *\  ,  t)  as    -T-*  0+ 
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CcIo(l/aXj  +  C,I,(2f5fJ 


(n-ar*...) 


l^aVt    e^nUb.-bO^' 


b,bitH  +  "    <fr* 


^ar 


-  isfi ... 


Thus 


|-  b.bi^-h^  (Zbi-bz.) 


+  b.  _2  +  ^« 

bi  *{ 


And  similarly 


73) 


-b,1 


r-*<r 


-M 


=  e 


Z  ,        -L 


xbl+lbfoKfl-cC  b2  ^  +  *  b2  b^ 


(74) 


Since  the  solutions  obtained  in  (71)  and  (72)  are  complicated,    they  are 
organized  so  that  the  character  of  the  solution  can  be  seen  as  clearly 
as  possible.     Equation  (70)  shows  that  there  is  no  disturbance  ahead 
of  the  wave  front.     It  is  also  noted  that  at  the  wave  front  (  "{     =  t-*|  =  0) 
0.   and  T.   are  singular,   both  containing  delta  functions  and     Q  Av\  ,  t) 

being  dominated  by  the  derivative  of  the  delta  function.      Thus    ©^(^  ,  t) 

2 

is  not  of  order  unity  and    £    ©j    is  not  small  in  comparison  to      Q    {  \  ,  t) 

at  the  wave  front. 
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Although  it  is  evident  that  the  labor  of  obtaining  higher  order 
corrections  (that  is,    higher  order  terms  in  the  perturbation  expansion) 
would  be  most  formidable,    it  should  be  equally  evident  that  such 
additional  corrections  would  have  a  nature  similar  to  that  of  the 
corrections  already  obtained;  namely,    they  would  exhibit  discontinuous 
behavior  at  the  wave  front  and  would  fail  completely  ahead  of  the  wave 
front  (that  is,    they  would  be  consistent  in  predicting  identically  zero 
temperature).     However,    it  is  to  be  expected  that  they  would  possess 
calculable  limiting  values,    from  the  left,   as  the  wave  front  is 
approached.     In  effect,   we  thus  have  what  appears   to  be  a  perfectly 
valid  way  of  determining  the  behavior  behind  the  wave  front,    but  one 
which  fails  at,    and  ahead  of  the  wave  front.     The  nature  of  the  failure 
of  our  original  perturbation  expansion  is  characteristic  of  singular 
perturbation  problems  [9],    and  we  are  therefore  led  to  seek  a  singular 
perturbation  expansion  which  will  resolve  the  difficulties  which  have 
been  encountered.     This  is  done  in  the  following  subsection. 

C.      SINGULAR  PERTURBATION  ANALYSIS 

We  have  seen  that  the  perturbation  solution  represented  by  the 
series  (12)  and  (13)  involve  no  difficulties,    except  those  of  analytical 
complexity,    for  any       %,  y  0       ,    that  is,    at  any  point  behind  the  wave 
front.     However,    this  solution  is  useless  for  points  ahead  of  the  front 

X  <o  ,  since  it  predicts  zero  temperature  for  fluid  and  solid  at  all 
points  ahead  of  the  wave  front;  physical  considerations    indicate  that 
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diffusion  through  the  wave  front  will  result  in  non-zero  temperatures 

ahead  of  the  front.     A  mathematical  indication  of  the  difficulty  may  be 

2 
seen  from  the  fact  that      £     q.        is  not  small  at  the  wave  front,    and 

4 
presumably  the  same  thing  could  be  shown  for      £    q       ,    etc.     Accord- 
ingly,  we  now  look  for  another  solution  which  will  be  valid  ahead  of  the 
wave  front.     To  do  this  it  is  convenient  to  introduce  the  new  variable 


Y  =    (n  -t)  /  e      =  x  -  Vt  (75) 

and  we  now  seek  new  solutions    0  (  V    ,  t),    T(  |  ,  t)  which  will  be  valid 
for    Y>o    • 

Making  the  change  of  independent  variable  from  (  1  ,t)  to    (Y  ,t) 
provides  the  differential  equations 

*  !*L£     «£§    +  b,0-T)  (76) 


*  c<    5j=~|I  +  <r  H-el^e-T)  (77) 


Note  that  the  transformation  (75)  changes  our  viewpoint  so  that,    in 
effect,   we  are  now  moving  to  the  right  with  the  front  so  that  the  fluid 
appears  to  be  stationary  and  the  pipe  appears  to  be  moving  to  the  left 
with  velocity  V. 

Note  that  in  contradistinction  with  the  situation  dealt  with  pre- 
viously,   now    £  (to  the  first  power)  appears  in  the  governing  equations, 
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2 

Accordingly    £   ,    not      g    ,    appears  as  the  natural  parameter  with  respect 

to  which  we  hope  to  construct  a  perturbation  expansion.     Thus  we  deal 
with 

(Later  when  we  "match"  these  new  solutions  to  the  solutions  discussed 
in  the  previous  subsection,    we  will  add  an  asterisk  to  distinguish  the 
new  solutions. 

The  conditions  that  all  these  new  functions  must  satisfy  are 

Si  — v   o  )      "Hi— *  o  as  ^-*oo  (80) 

0/sTt'so  for  t  =■  o  (81) 

for  i  =  0,  1,  2,    ... 

Additionally,    they  must  satisfy  appropriate  conditions  at  the  wave 
front.     Actual  selection,    or  determinations,    of  these  conditions  for 

£   =    0  will  be  made  in  the  next  subsection  where  the  matching  process 
will  be  treated. 

Substituting  (78)  and  (79)  into  (76)  and  (77)  comparing  powers  of   £    , 
gives  rise  to  the  following  sequence  of  equations.      For  the  fluid,   we 
will  have 
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*   9*    =    e0  +  b,  (0,-To)  (82) 


o^e*  =  e,  +  b,  (e,  -I) 
*  el    =  e,+  b,  (&2-7i) 


(83) 


(84) 


and  for  the  solid  we  obtain 


/ 


L     =0  (85) 


*mT7  =-l,  +•  l«-b2  (ea-Toj  (86) 

^Ti"  =-TL  +7i  -  bi  (0,-TI)  (87) 

it' 


•      •      • 


(where  the  prime  denotes  partial  derivative  with  respect  to      ^      and 
dot  denotes  partial  derivative  with  respect  to  time). 

Condition  (80)  and  equation  (85)  show  that 

To~    0  (88) 

Now  consider  the  first  non-trivial  equation  in  the  sequence, 
namely  (82).     This  is  a  classical  diffusion  equation,    the  solution  of 
which,    subject  to  (80),    is 

-b,t 
*.f?'*)  =Ao<?     '  er*c  (1/2 /^?  J  (89) 
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The  constant  AQ  is  available  to  satisfy  a  boundary  condition  for  ^  =0. 
The  next  equation  to  be  considered  is  (86)  which  in  view  of  (88)  may  be 
written  as 


Tt  =  -  b2  ec 


(90) 


Together  with  conditions  (80)  and  (81),    we  may  write  down  the  solution 
of  (90)  as 

bt    °° 
T\$,t)  =  A0  b2  e  '     Terfc  (?/2{S£t  )  ^  (91) 

f 

Note  that  no  new  constant  of  integration  is  required  in  the  last  step. 

The  next  equation  is  (83)  which  is  similar  to  (82)  except  that  T    , 
given  by  (91),    is  involved.      The  solution  could  be  obtained  but  we  will 
not  do  more  than  remark  that  another  constant  of  integration  will  be 
introduced.      Then  (87)  provides   T£;  no  new  constant  is  introduced. 
Next  (84)  provides  a  means  of  determining     Q  2'>  another  constant  of 
integration  is  involved. 

Briefly  the  solution  for     Q  •   (i  =  0,  1,  2,  ...  )  involves  a  constant 
A-,   which  is  actually  a  function  of  time  t.     Since  the  equations  for  T. 
are  first  order,    their  solutions  do  not  require  the  introduction  of 
arbitrary  constants  or  functions,    equation  (80)  and  (81)  providing  the 
necessary  conditions. 

No  difficulties  are  encountered  (except  those  of  increasing 
analytical  complexity)  in  continuing  this  process  and  we  expect 
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convergent  series  (78)  and  (79).     However,    as  yet  we  have  made  no 

move     to  relate  these  expansions  to  any  conditions  at  the  left  which 

might  serve  for  the  establishment  of  the     as  yet  undetermined     functions 

A  (t).     We  proceed  to  do  this  in  what  follows  immediately. 
i 

D.      MATCHING  OF  SOLUTIONS 

We  now  employ  an  asterisk  to  distinguish  the  solutions  (  0  .    and  T.   ) 
of  the  immediately  preceding  subsection,   valid  ahead  of  the  front 
(i.  e.  ,    for     «  >  o   ),    from  the  solutions  (  Q  ■  and  T. )  of  the  earlier 
subsection,   valid  behind  the  front  (i.  e.  ,    for     X.  >*>).     We  will  employ 
the  following  notations. 


(92) 


•&■    Ti(*,t)   *Tc  it)  =  Ti  (93) 


[*V>]      - 


=       e?(t)     J/  (94) 


?*0 


Tito>*)\         =      Tl    (*)      =Tc  (95) 


(Note  that  here  an  overbar  does  not  indicate  a  Laplace  transform.  ) 
Examples  are 

e0    =  e"b,t  (96) 
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TQ     =    0 


(97) 


=     0 


=  ka  e 


-b,t 


CD 


(98) 

(99) 
(100) 

(101) 
(102) 


The  functions  Q-     and  T  .     contain  as  yet  undetermined  functions  of  t, 
namely  A     =  A    (t),   A,    =  Aj(t),    etc.  ,    and  we  must  be  able  to  evaluate 
these  systematically.     Our  procedure  is  to  evaluate  so  that  the  per- 
turbation sequences  to  the  left  and  right  of  the  front,    truncated  so  as 
to  include  only  terms  containing    €  to  an  arbitrary  even  power  (and 
lower  powers)  will  match  at  the  interface. 
The  first  matching  is  to  require 


-  -  * 

a  -  *> 


(103) 


To  =t;* 


(104; 


Equation  (104)  is  automatically  satisfied.     Equation  (103)  provides  a 

condition  which  can  be  used  to  determine  the  constant  A    .     We  find 

o 
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A     =  A    (t)  =   1.      Thus,    our  first  (matched)  approximation  is  given  by 
(29),    (31),    (89),    and  (91)  with  Aq  =   1  in  the  latter  two. 
The  second  matching  is 


0o  +■  €*"  0,    =    Bo   +  i  9  *  +   fe*  a2 


(105) 


To  +  6  T,    =  T0    +€lT     +    6   T2  (106) 

*  *         *  * 

The  functions      9|   ,     9-  ,    Ti  ,    and  T?    have  introduced  two  additional 

unknown  functions  A^(t)  and  A    (t).      Equations  (105)  and  (106)  provide 

conditions  for  their  determination. 

Similarly  the  next  matching  is 


-  2.  -  4_  _£  _1f.        2-rt  3  -#        t   _  * 

0O+  6   0,  +  6  9t   =    90    -t-  £  a,   +  £  £z  -+  €   e3   +  £   &«  (107) 


-  2.-  <f-  -*  -*         2-*  3-*  *— * 

To-*-^77+^T2=X   +    6T;-h£Tz     h-£T3    -+  £  T^.  (108) 


in  which  two  additional  functions,    A-^(t)  and  A.(t)  appear  and  which 
afford  means  for  their  determination.     Proceeding  in  this  manner,    all 
functions  A.(t)  may  be  determined.     At  each  stage  we  have  matched 
solutions  to  the  same  degree  in    £   ,    on  either  side  of  the  front. 

By  this  procedure,    eventually  (1)  and  (4)  are  satisfied.      Moreover, 
at  any  stage,    the  necessary  coefficients  A.(t)  are  calculable  so  that 
the  approximation  affordytruncated  series  ,     available  both  to  the  left 
and  to  the  right  of  the  front,    with  matching  at  the  front.      Thus  it  appears 
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that  we  have  succeeded  in  describing  a  process  which  provides  two 
sets  of  series  expansions,    one  valid  behind  the  front  and  the  other 
valid  ahead  of  the  front.     (A  "    et"  consists  of  a  solution  0    for  the 
fluid  and  a  solution  T  for  the  solid.  )     The  labor  of  constructing  the 
successive  functions  that  occur  in  these  series  is  considerable  and 
we  have  actually  obtained  only  a  few  of  them.     We  have  actually  made 
the  "first  matching"  and  could,    with  additional  labor,    determine 
sufficient  additional  functions  to  accomplish  second,    and  even  third 
matchings.     However,    the  writer  confesses  to  a  feeling  of  uncertainty 
in  this  entire  matter  of  constructing  the  ordinary  and  the  singular 
perturbation  expansions  and  the  matter  of  performing  a  matching 
where  it  is  necessary  to  join  them.     Accordingly  he  suggests  that  it 
may  be  fruitful  to  consider  this  matter  further,    employing  the  insights 
which  might  be  provided  by  a  more  sophisticated  mathematical 
treatment. 

E.      A  DIFFERENT  SOLUTION  AHEAD  OF  TEMPERATURE  FRONT 

The  solution  given  in  this   subsection  may  be  of  some  interest 

i 

although  it  involves  a  certain  basic  inconsistency.     In  this  solution 
we  are  concerned  with  the  temperature  variation  ahead  of  the  front 
and  are  able  to  consider  both  radial  and  axial  variation  in  the  fluid, 
as  well  as  both  axial  and  radial  diffusion  in  the  fluid.     However,   we 
assume  that  the  temperature  of  the  solid  is  uniformly  zero  and  that 
there  is  perfect  heat  transfer  through  the  interface;  that  is,    we 
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consider  an  isothermal  boundary  condition  for  the  fluid.      (It  would 
not  alter  the  problem  profoundly  if  a  finite  surface  coefficient  of  heat 
transfer  were  assumed.  )     Thus,    this  problem  is  not  a  conjugated 
problem.      For  a  left  end  boundary  condition  we  assume  that  the  fluid 
temperature  is  a  constant,    independent  of  radius,   which  depends  on 
the  location  of  the  front.     Accordingly,    the  degree  of  approximation 
and  inconsistency  in  this  solution  is  considerably  greater  than  is  to 
be  found  in  other  problems  treated  herein.     However,    the  results  may 
be  useful  in  considering  the  temperature  distribution  ahead  of  the 
front  at  quite  small  values  of  time,   before  significant  radial  variation 
has  built  up  in  the  front. 

We  will  consider  equation  (II- 1)  and  (II-2)  with  isotropic  thermal  con- 
ductivities   and  subject  them  to  a  change  of  variable 

§«x-Vt  ;  t-t  (109) 

Then  equation  (II- 1)  reduces  to 


ti  +L~bi+'bJ?    -  _L  12  (no) 


Similarly  equation  (II-2)  reduces  to 

ti  + 1  £T  +  ti  +  v  *T  =  i  M" 

We  note  that  the  effect  of  the  change  of  variable,    (109),    is  to  change 
our  view  point  as  follows;  the  fluid  remains   stationary  and  the  pipe 
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wall  moves  to  the  left  with  the  same  velocity    V    of  the  fluid.     If  the 
velocity    V    is  large  then  we  state  the  following  assumption 

Assume  that  for     V  >  o     the  pipe  wall  is  essentially  at  zero 
temperature. 

Hence  we  will  concentrate  our  attention  on  equation  (109)  and 
find  the  solution  when  equation  (109)  is  subjected  to  the  following 
conditions 

0=o  ^-+oo  ,  o<r*Ct  ,  t>o  (113) 

9  =  o  0^£<.co  >  fa  OL  ,  *>o  (114) 

^=o  o<-t  <co  j  r=o  ,  -t>o  (115) 

9=o  o<^<a>  j  o  <r  <.  a  ,  t-a  (116) 

(Note  that  r=0     implies      x=\/t,    and      ^  — *-  oo        implies  x  — »  °°     ) 

Assume  the  solution  to  be 

oo 

9   =   Ykn  tti'VJA^r)  (in) 

condition  (114)  yields 

Jo(Xnaj  =o  (118) 


and  condition  (115)  is  already  satisfied. 
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Next,    taking  Laplace  transforms  of  (110)  and  (117),    we  find 

tin  -(X.+S-)  t  =o 

A  solution  of  (119)  which  satisfies  (113)  is 

Thus,   writing    Cw=sAnBJl    ,   we  have  * 

00  r 


»;/ 


Condition  (112)  then  yields 


oo 


i-    £CJ.M 


Using  the  orthogonality  of  Bessel  functions,    viz, 


a 


/jtMJi(XIBr)rJr*-^^  (XMaj 


if  n  =  m  and  zero  if  n  ^  m,    we  find, 


r    -     '      z 


5      W       Jt(W) 


Hence 


(119) 


(120) 


(121) 


(122) 


(123) 


(124) 


oo 


n=i 


(125 
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Inverting  (125)  we  find 


oo 


^fW^  ir* 


(Note  ;  that  solution  (126)  is  for  the  case  when  the  wave  front  is 
maintained  at  unity.  )    It  is  now  possible  to  use  Duhamel's  super- 
position integral  when  the  condition  at    "P  =  0  is 


-b,i 


(127) 


which  is  the  simplest  estimate  we  have,    Cf.    (96).     Thus,   we  obtain 

t 


8(l>0*j.     /C        i 


2»t-S 


(128: 


where  the  term  in  brackets  may  be  obtained  by  differentiating  (126). 

Obviously,    this  solution  is  for  cylindrical  geometry,    but  the  case 
of  slab  geometry  differs  only  in  detail. 
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V.     CASE  WHERE  FLUID  TEMPERATURE  IS  AVERAGED  AND 
AXIAL  CONDUCTION  IN  CONDUIT  IS  NEGLECTED 

A.     SLAB  GEOMETRY 

1.       Development  of  Solution 

In  this  section  equations   (11-19)  and  (11-20),    together  with  the 
appropriate  boundary  and  initial  conditions  will  be  further  developed. 
For  convenience  of  reference  these  equations  are  rewritten  below. 


bZa     V  ~69    _    I    "6&    ,      « 


L*-V.  !£.   =i-  1£  +.  2L    (q.j)  (1) 

bx2-  <*  ax     <*   dt     a»cf  *       ' 


&      SiS  (2) 

*3*     «*  at 


0  =  o                       oix<a.  ,  t  =o  (3) 

0*1                             X  =  o  ,  t>o  (4) 

0->O                           X-*<»  >  t^O  (5) 

J  =  C>                      o<£  <  b  ,  t  =  o  (6) 

2"  =  o                          y  =  b  ,  *£o  (7) 

-K|M(e_TJ       y=<?  j  t>fl  (8) 

Note,    that    9=9    (x,  t)  appearing  here  is  the  average  fluid  temperature 
as  defined  in  (11-18);  we  have  omitted  the  asterisk  uniformly  throughout 
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this  section.     Also,    T  =  T  (x,  y,  t)  is  the  temperature  in  the  duct  wall; 
we  have  omitted  the  subscript  2  on  y;  in  this  section,    the  symbol    y 
denotes     y      in  figure   lb.     We  also  note  that  this  analysis  neglects 
the  axial  conduction  in  the  duct  wall.      Finally,    note  that  the  surface 
heat  transfer  coefficient  appearing  above  is  the  "equivalent"  film 
resistance  coefficient   £    discussed  following  (11-17);  we  omit  writing 
the  tilde  in  (1)  and  (8)  above  and  in  what  follows. 

In  order  to  solve  (1)  and  (2),    Laplace  transformation  is  used. 
Equations  (1),    (3),    (4)  and  (5)  become 


•-  ffij,.. 


(9) 


^5"  X=0  (10) 


Q-^-o  X—*oo  (11) 


Similarly  equations   (2),    (6),    (7)  and  (8)  become 

^f  *-i-T  (12 


-**^»-K(*-t)  o<x<co       ,       ^o  (13) 

<£T  =  o  o<*{oo      ,       y^b  (14) 

^  * 
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The  solution  of  (12)  is 


T=  Ae   +Bed 


(15) 


where 


y-ik 


(16) 


Condition  (14)  permits  evaluating    B    in  terms  of  A,    and  we  arrive 


at 


f  -  A 


f  W     i(*fc-a)l 

e  +  e 


(17) 


Since  the  terms  in  parenthesis  are  clearly  independent  of  x,    the 
imposition  of  condition  (13)  implies  that  A  in  (17)  must  be  a  function 
of  x  and  s.     Thus  using  (13)  to  express  Q  in  terms  of  T,    we  find 


9  =  pk).  A(x7s) 


(18) 


where 


/»«).-£  j(.-2>(l  +  !*b) 


(19) 


Substituting  from  (17)  and  (18)  into  (9)  yields  a  homogeneous 
second  order  ordinary  differential  equation  for  A(x,  s).     After  re- 
arranging and  substitution  for      /3  (s)  this  becomes 


(20 
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where 


3(s) 


K, 


) 


(2i; 


The  general  solution  of  (20)  is 
A(x,s)=Ce*p 


i  r  c  f 


122.) 


where  C  and  D  are  constants  of  integration  which  depend  on  S.     The 
boundary  conditions  (10)  and  (11)  in  conjunction  with  (18)  define  the 
following  conditions  for    A  (x  ,s) 


A  (x,s)  >  0 


x  =  0 


oo 


(23) 
(24) 


Condition  (24)  clearly  implies  that  D  S  0  whereas  (25)  yields 


c    = 


s./»(0 


(25) 


Therefore 


A(x,s)=  -!— ■  ^p 


2^      f 


(26) 


Substitution  of  A(x,  s)  from  (26)  above  into  (18)  and  (17)  yields  the 

fluid  and  pipe  wall  temperature  distribution  in  Laplace  transform  form; 


viz. 
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0(x^)= 


i..e*p 


$  W(b^) 


(27) 


and 


f(x,j|,s)«< 


c^h(b-a)] 


b  n 


(28) 


T^ 


MW(bf) 


=  f(^).a(x,s) 


In  order  to  invert  (27)  and  (28),    we  first  consider    Q(x,  s).     Denoting 
the  radicand  by  m     we  have. 


m 


.  {Mi i : w w  w fejj&fridk 


i+&(kl)W(H) 


(29) 


u 


It  will  be  convenient  to  cast  this  equation  in  another  form.     Write 


x=  -i 


)  - 


"I 


m 


(30) 


so  that 


S    =  -*  x 


«\ 


I 


=  a 


(31) 


Then  (29)  becomes 


m  = 


„  ^kkf-^ty^AttS'+h-izA 


B„  _bX  to-  (bX; 


Ml 


(32) 
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where 

Bm=  T-  ;        V  "k7  (33) 

are  Biot  numbers  for  the  wall  and  for  the  fluid,    respectively.     We 
are  concerned  with  finding  the  singularities  of  (27).     Except  for  the 
pole  at  s=0,    these  will  all  be  branch  points  corresponding  to  the  zeros 
and  poles  of  m   .      Considering  the  numerator  and  denominator  of  (32) 
separately,   we  first  note  that  the  denominator  has  poles  for       da    = 
(n+  /2)tf     ,  n  =  0,  1,  2.  .  .     .     These  will  also  be  poles  for  the 

numerator  unless 

We  may  assume  that  this  equality  does  not  hold  for  any  integer    n 
for  otherwise  we  could  make  an  infinite simally  small  change  in  any 
of  the  physical  parameters  so  as  to  assure  that  the  equality  is  violated, 
without  making  any  significant  change  in  the  physical  problem.      Thus, 
the  only  zeros  and  poles  of  m     correspond  to  the  zeros  of  the  num- 
erator and  denominator,    respectively. 

The  condition  for  zeros  of  the  numerator  (i.  e.  ,    zeros  of  m    ),    is 
that 

tan(Xb)      =f.    (a)  (35) 


where 
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«(*)- 


(V.    \l     %2>  X2 

B*       \ioCf)       <*+  ^36) 


0*       ^ 

and  the  condition  for  zeros  of  the  denominator  (i.  e.  ,   poles  of  m   ), 

is  that 

tan(bX  )  =f2  (   X  )  (37! 


where 


£,(   A  )      =     JLb»-  (38) 


bX 

The  real  zeros  of  (35)  and  (37)  may  be  found  easily.     Their  location 
is  illustrated  in  Figure  3.      The  relative  position  of  the  dashed  line 
and  the  dot-dash  line  in  Figure  3,    is  important.      For  small    X   it  is 
evident  that  f   (    X     )       <(    f2(     X    )•     Also  f   (    X    )  has  a  vertical 
asymptote  for  a  value  of     X    greater  than  that  of  its  single  zero.      Thus 
the  positions  are  as  shown  in  Figure  3. 

Finally,    we  must  assume  that  there  are  no  other  zeros  than  the 
real  zeros  indicated  in  Figure  3.     However,    both  (36)  and  (38)  satisfy 
the  hypotheses  of  Appendix  B  which  shows  that  there  are  no  non-real 
zeros. 

Thus  it  is  clear  that  the  only  singularities  of  m     are  poles; 
furthermore,    m     has  only  real  zeros  and  poles,    infinite  in  number, 
and  alternating  with  each  other,    of  which  that  having  the  smallest 
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o  - 


y   =   tan(bX) 

y   =    fx(\) 

y  =  Vx) 

zeros,  x  -  poles 


Figure   3:      Poles   and   Zeros  of     m(s) 
(Slab   Geometry) 


118 


numerical  value  is  a  zero.     If  we  designate  the  abcissae  of  the  inter- 
sections denoted  by  0  in  Figure  3  by  the    symbol    \      and  those 
denoted  by    x    by  the  symbol       Xn  »   we  conclude  that  the  singularities 
of  (27)  are: 

a.  branch  points  (corresponding  to  zeros  of  m   )  at 

S  =  ~~  o(      \  ,    and 

m  rt  n        ' 

b.  branch  points  (corresponding  to  poles  of  m    )  at 

S    =    '     *  An2 

m    /     n 

c.  a  pole  at  S  =  0  with  residue     1 

Furthermore,    the  branch  points  of  type  a  and  type  b  alternate  with 
each  other,    that  closest  to  the  origin  being  of  type  a. 

Thus  the  hypotheses  of  the  inversion  theorem  given  in  Appendix  A 
hereof  are  satisfied.     Accordingly,    we  may  immediately  write  the 
inverse  as  follows 

*k 

It 

<K  (39) 

i        L  '       J 


00        it 


0(x,t)  =   1   -    -Le^(^)V    /    £_   ^    x//(l) 


where       ^Cl)   =     (l-*w(-l)|  ) 


Similarly  T(x,  y,  s)  as  given  by  (28)  is  conveniently  organized  in  the 
form  of  the  product  of  two  Laplace  transforms.     Accordingly,    its 
inverse  may  be  easily  obtained  in  the  form  of  a  convolution  integral. 
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oo 

r — ■  .     »2. 


(40) 


where         X«        are    the  roots  of  (37)  and  the  integrand  contains  the 
function  (39).     We  have  used  the  formula 


-/ 


Cost  (OC^) 


CaSt^'i)  +4{bSl'i)  /^(bs'/ij 

XM    <f  n  Ax«(\na)Cos(\n\>) 


(41) 


Z-.     2XMb     +  /U^  (2bXn) 


which  may  be  derived  directly  or  inferred  from  equations  given  in 
[2],    examples  5-3  and  7-11. 

2.       Asymptotic  Approximations 

Thus  the  desired  solution  of  the  system  described  by  equation 
(1)  through  (8)  has  been  achieved.     We  will  not  in   this  work  discuss 
computational  aspects,    but  will  now  obtain  some  asymptotic 
app  roximations . 

Equation  (27)  can  also  be  written 


0(*,s)^le*p 


%-*MH^-^y*i---) 


(42) 
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where 


X^W(fcj)  =  5A_ 


sf    sV 


3,6 

2.  S  t 


n  s  6 


+  —  -rrr  -  —  -7i  + 


3«i      '5    *>       35    ** 


m 


(43) 


Consider  the  radicand  to  be  expanded  in  ascending  powers  of  s. 


e(x^)  =  4-«*p 


2* 


f 


1     *n 


(44) 


Considering  small  values  of  s,    which  corresponds  to  large  t, 
the  first  asymptotic  approximation  is  obtained  by  neglecting  all  powers 
of  s  except  the  zeroth  power.      Correspondingly 


0(x,s)    =    i 


(45) 


which  yields 

eu,t;  =  i 


(46) 


This  is  evidently  correct  as  an  asymptotic  expansion  since  for  any 
finite  x  the  temperature  approaches  unity  for  sufficiently  large  time, 

A  nontrivial  asymptotic  expansion  is  obtained  by  retaining 
the  first  power  of  s.      Correspondingly 


a(x,s)^  i  **p 


(47) 


The  reason  for  writing  Am  /Af  ,    the  area  ratio,    in  place  of 
b/a  which  naturally  would  appear  in  the  last  term  of  the  radicand  will 
appear  in  Section  VB2. 
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The  inversion  of  (47)  is  tabulated  [29],    thus 


0(*,t)«     1 


^  0$)  ♦ e"  ■*  (ffi) 


(48) 


where 


v*-*/(i+&c*AJf;-%+t) 


(49) 


is  a  retarded  velocity  which  is  identical  to  that  obtained  by  Munk  [24], 
and 

*         I  c 1 ' 

(50 


1*=  l//^H/fS 


is  an  effective  thermal  diffusivity  for  the  fluid.      This  solution  should 
be  compared  with  that  given  in  (VII- 10)  where  a  different  set  of 
idealizations  leads  to  a  similar  result.     There,    the  retarded  velocity 
(VII-5)  is  the  same  as  (49),    but  the  diffusivity  (VII-6)  is  slightly 
different  from  that  given  in  (50). 

For  the  corresponding  solution  for  the  conduit  wall  temp- 
erature we  employ  (48)  in  (40). 

Retaining  the  second  power  of  s  would  yield  a  useful 
improved  asymptotic  approximation  if  one  could  find  a  simple  explicit 
inversion  formula.     Unfortunately  nothing  simpler  than  the  general 
inversion  given  in  Appendix  A  has   come  to  light,    and  this  has  provided 
the  complete  solution  given  in  (39)  and  (40). 
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B.      CYLINDRICAL  GEOMETRY 
1.       Development  of  Solution 

A  similar  treatment  to  that  given  in  Section  VA  is   employed. 
Only  the  significantly  different  details  will  be  treated  explicitly. 
Equation  (2)  will  be  replaced  by 


?LT  +.  1ST  =  ±  "£T 


(51) 


with  the  corresponding  conditions 

T=  o  '     O.  £f <b  )        t=a  (52) 

21  =o  f=*>  J         *>°  (53) 


M<>r        K        J 


(54) 


replacing  (6),    (7),    and  (8).     Remarks  similar  to  those  following  (8) 
are  applicable  here. 

As  previously,    Laplace  transformations  are  made  with 
respect  to  t.     The  solution  of  the  transformed  version  of  (51)  is 

f  =    Al0(|rJ+8K^r)  (55) 

where 


j  =  {VZ 


(56) 


is  the  same  as  previously.     Here  I     and  K     denote  modified  Bessel 
r  J  o  o 

functions. 
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Proceeding  as  in  VA1,    one  arrives  at  equations  similar  to 
(27)  and  (28)  except  that  the  function  tan{|  (qb)  appearing  therein  is 
to  be  replaced  by  another  function  which  we  denote 

W(x,P)«      I.M*W-WK.(X)  (57) 

I.(px;K,(x)+K.(px)Ii(K) 

(The  initial  b  denotes   "Bessel.  "     The  remainder  (anh)  is  intended  to 
indicate  a  similarity  to  the  function  tan£(x)      .  ) 

Corresponding  to  (29),    we  have  in  this  case,    the  equation. 


_  (&jUt+(^W4)fe(^  frjfe] 


m  =  v-^£ 


Corresponding  to  the  well  known  relation 


(58) 


(59) 


it  is  not  difficult  to  show  that 


tkM£(Xjfi)     =i     bfl^i  (I*,  p)  (60) 


where 


(61) 
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Thus  corresponding  to  (32),    upon  defining 

*=-^  (62) 

we  obtain 

w.  -  Ew  [(£)  -  fol'N  MbMA)  fe  ^g]  (63) 

8W   -(bx)  bom(U,a/b) 


where 


d   ...  b*  .  R*-  Z^ 


(64) 


(Note  the  factor  2  in  the  present  definition  of  B/    ) 

From  this  point,    the  development  is  precisely  the  same  as  in  VA1. 
Considering  the  counterparts  of  equations  (35)  and  (36),    precisely  the 
same  treatment  is  applicable.      Figure  4,    shows  a  graph  similar  to 
Figure  3.      Figure  4  was  actually  constructed  by  computer  analysis  in 
order  to  show  that  the  function  ban  (  X  ,  fi  )  behaves  like  tan(x)  ;  the 
value  of    ft    used  in  constructing  Figure  4  was     (I  =    U.9.  . 

There  is  one  difference  which  we  would  prefer  to  ignore. 
In  the  previous  development,    the  fact  that  m     had  no  non-real  zeros 
or  poles  was  proved  by  the  argument  contained  in  Appendix  B.      To 
construct  a  similar  proof  for  the  case  in  which  tan  x  is   replaced  by 
ban  (x,  (2  )  appears  terribly  formidable  and  we  simply  omit  doing  so. 

Thus,   a  solution  corresponding  to  (39)  and  expressed  in 
exactly  the  same  notation  as   (39)  (that  is,    identical  to  (39))  is 
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o  - 


y  =   ban(bX,3> 


y  =    f,(X) 

y  =    f2(X) 

zeros,   x  -  poles 


Figure  4:     Poles  and  Zeros  of     m(s) 
(Cylindrical    Geometry) 
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applicable  in  the  present  case.     However,    the  function  F  (r,s)  in  (28) 
must  now  be  replaced  by 


R<-,5)  = 


(65) 


The  inversion  of  (65)  may  be  obtained  using  a  solution  discussed 
in  Section  13.  4  of  [7],   by  first  specializing  to  the  case  of  a  unit  step 
input,   and  then  by  differentiating  the  result  with  respect  to  time.     The 
result,    corresponding  to  (40),    is 


where    yuw     are  the  roots  of 
and 


(66) 


(67) 


;£M-  Jk(/,x^Yw(P8J-Yk(NJm((98) 


(68) 


2.       Asymptotic  Approximations 

It  is  possible  to  proceed  here  in  exactly  the  same  way  as  in 
Section  VA2.  The  principal  difference  is  that  in  the  present  case  of 
cylindrical  geometry,  in  place  of  (43)  we  now  define  ~X-  differently 
and  arrive  at  different  coefficients  in  the  expansion.      Thus 
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X.-(bt)b*.B(H,«/b)  =  ~a^-<^)  *+■•■  (69) 

where  .  .  .    denotes  terms  in  s    ,    etc.     It  requires  painstaking  work 
to  establish  that  the  coefficient  of  s     vanishes  and  that  the  first 
power  term  is  as  in  (69)  but  the  calculations  are  not  illuminating  and 
are  not  given  here. 

As  in  the  case  for  slab  geometry,    the  first  non-trivial 
asymptotic  expansion  is  obtained  by  retaining  the  first  power  of    s 
in  (69);  correspondingly,   we  arrive  at  equation  (41)  the  inverse  of 
which  is  given  in  (48).     In  writing  (47)  we  anticipated  the  present 
development  by  writing  A^^  /A £   in  the  final  term  of  the  radicand  so  as 
to  be  able  to  accommodate  to  both  slab  geometry  and  cylindrical 
geometry  with  one  equation. 
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VI.     A  GENERALIZATION  OF  A  PROBLEM  CONSIDERED 

BY  ARPACI 


This  case  is  described  by  (11-10)  and  (11-14)  subject  to  conditions 
that  will  be  stated  below.     The  problem  description  is  thus 


(1) 


?t=fh  (e~T)  (2) 

6  =  1  X=«o  j  t£o  (3) 

6  =  0  X-^  co  ^  t>o  (4) 

0  =0  °<X  <oo  j  t  -o  (5) 

T-  0  0  <  X  <&  j  *~o  (6) 

0-0  (x,  t),    T  =  T(x,  t)  are  average  temperatures  in  the  sense  of 
(II-7),    (II- 8),   but  we  omit  using  the  asterisk  which  distinguished  this 
significance.     The  heat  transfer  coefficient  h  in  (1)  and  (2)  above  is 
defined  by  (II-9). 

This  problem  is  a  generalization  of  that  considered  in  Section 
7-5  of  [2],    the  generalization  coming  about  by  the  inclusion  here  of 
axial  conductivity  in  the  fluid,    a  mechanism  not  considered  in  Arpaci's 
analysis.     This  problem  is  a  specialization  of  the  problem  considered 
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in  Section  IV  hereof  since  here  we  neglect  axial  conduction  in  the 
duct  wall,   whereas  this  mechanism  is  included  in  the  analysis  given 
in  Section  IV. 

Taking  Laplace  transformation  with  respect  to  time  t,    one  obtains 


ST  =  bi(e-T) 


where 


b,  = 


it 


b». 


-(IP 


Combining  (7)  and  (8),    one  obtains 


A- 


V 0      _V_  "be    _  /  £.  ,  L         t>|ba.   \ 

'  ^  bx    uf        s+b2; 


W 


e  =o 


subject  to  the  conditions 


5-  1 


x  =0 


(7) 

(8) 


(9) 


(10) 


(11) 


0 


oo 


(12) 


Thus,    the  solution  of  (10)  may  be  written 


0=    1  exp 


%-W*V*- 


bik 

S+bz 


(13) 


an  increasing  exponential  having  been  rejected  in  order  to  satisfy 
(12).     We  write  (13)  as 
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e  ^ 


S         r 


(14) 


where 


x=f( 


— 7 » 


0 


(15 


f-2  0- 


) 


(16) 


CO 


t      4*f 


(17) 


It  is  easy  to  show  that     A  >  i-O   0        •     However,   we  want  to  indicate 
the  magnitudes  of     X     and      U      relative  to  b~  .     To  see  this,    consider 
the  function 


f">"5;-<V,-ty*SS 


(18! 


the  vanishing  of  which  corresponds  to  the  roots  p;  \       ;  viz. 

f(-  X     )  =  f(-  p     )  =  0.     We  note  that  the  coefficient  of  the  quadratic 

term  is  positive,    and  we  note  that 

f(-b2)  =  -bXb2      <     0.  (19) 

Thus,   we  easily  see  that, 

X>t>i>p  >   O  (20) 

Thus  the  condition  of  the  theorem  stated  in  Appendix  A  are  satisfied, 
and  we  can  write  the  inverse  of  (14)  as 
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where 


eix,*)-ir£*p($L)- 


bi 


i* 


^(xfwpn     (2i) 


00 

r  -a-fc 


'<*.)- ^g^ 


(22) 


Using  the  convolution  theorem  to  invert  (8),   we  obtain 


r 


(23) 


We  should  remark  that  equations  (1)  and  (2)  above  may  be 
obtained  by  specializing  equations  (61)  and  (62)  of  Hiep's     treatment 
[13].     However,   we  note  again  that  Hiep's  results  depend  upon 
numerical  inversion  of  Laplace  transformations,    a  process  which 
our  experience  indicates  is  of  doubtful  accuracy  and  reliability. 
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VII.     AN  IMPROVEMENT  ON  A  RESULT  BY  W.    MUNK 

Under  the  assumption  that  radial  conduction  is  so  efficient  that 
essentially  the  same  temperature  prevails  at  all  points  of  a  particular 
fluid  and  metal  cross  section,    W.    Munk  [24]  obtained  a  very  simple 
approximate  solution  for  what  he  called  "The  Delayed  Hot  Water 
Problem.  "    In  the  present  section  we  obtain  a  very  simple  improve- 
ment over  Munk's  result.     Hiep  [13]  has  also  obtained  the  improvement 
given  here  and  he  mentions  that  it  has  been  obtained  by  others.     Our 
reason  for  presenting  it  is  that  it  is  a  natural  deduction  to  be  drawn 
from  developments  already  made  in  this  thesis. 

We  repeat  equation  (11-13)  and  the  corresponding  conditions. 

n,+  k,  Af  i  7>x*     V"bx     V+  kfA,  «J  at  (1) 

T  =  1  at  x  =  0        i       t  ^  0  (2) 

T  =  0  as  x— ►  «>      j        t  =  0  (3) 

T  =  0  at      O^X  <<x>       >        t  =  0  (4) 

We  have  used  T  to  denote  the  common  value  of  average  temperature. 

By  defining 


Jim  "m    *• 


--o^W^/o+f^) 


(6) 


133 


we  may  write 


Taking  Laplace  transforms  with  respect  to  time,    we  obtain 


fcx2         bx 


The  solution  of  (8)  satisfying  all  conditions  is 


T=  j  **? 


V*  IN* fT  S 


Inverting  (9)  by  the  use  of  standard  tables  [29]  we  obtain 


TW)=! 


fi*t  ) 


(7) 


(8) 


(9) 


(io; 


This  is  the  generalization  described  in  the  title  of  this  section. 
Munk's   result  is 

T     (x,t)  =     H       (x-V*t)  (11) 

where    H     (  )  denotes  the  Heaviside  unit  step  function.      Thus 

the  present  result  (10)  is  like  Munk's  in  predicting  a  temperature 
wave  front  moving  with  a  reduced  velocity  V    ,    but,    unlike  Munk's 
it  also  provides  for  a  softening  of  the  wave  front  due  to  axial 
diffusion  in  fluid  and  pipe. 
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VIII.     DISCUSSION 

A.     SCOPE  OF  ANALYSIS 

The  general  problem  as  formulated  by  equations  (II- 1)  and  (II-2) 
for  the  case  of  cylindrical  geometry  and  equations  (II-3)  and  (II-4) 
for  the  case  of  slab  geometry,    even  with  the  simplifying  assumptions 
of  subsection    IIA2,     are  mathematically  very  involved.     These 
equations  contain  two  dependent  variables,    namely  the  fluid  and  the 
pipe  wall  temperature  each  of  which  is  a  function  of  three  independent 
variables.     There  is  also  a  formidable  number  of  boundary  and  initial 
conditions  to  be  satisfied.     In  Section  III  the  problem  is   reduced  to 
the  solution  of  an  integral  equation  the  only  unknown  in  which  is  the 
temperature  at  the  interface,   which  is  a  function  of  only  two  variables; 
no  auxiliary  conditions  are  required.     However,    as  may  clearly  be 
observed  by  the  reader,    the  integral  equation  of  interest  is  math- 
ematically very  involved,    and  retains  a  great  share  of  the  complexity 
of  the  original  problem,    for  the  following  reasons. 

Firstly,    the  unknown  function  appears  only  under  the  integral 
sign  (the  so  called  integral  equation  of  the  first  type).     Secondly, 
there  are  two  integral  signs  under  which  the  function  of  interest 
appears.     Thirdly,    one  of  the  integrals  has  a  semi-infinite  range, 
namely  o      to      co       .      Fourthly,    the  kernels  of  integrations  are 
complicated  functions. 
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The  complexity  of  the  general  problem,    either  in  its  original 
form  or  as  compressed  into  this  integral  equation  is  such  as  to  have 
led  us  to  focus  attention  upon  further  simplifications  which  are  more 
tractable  from  a  mathematical  viewpoint  and  the  analyses  of  which 
are  treated  in  other  sections  of  this  thesis.     However,    the  integral 
equation  formulation  is  believed  to  be  of  significance  and  its  treat- 
ment,   probably  by  numerical  procedure,    is  considered  feasible. 

Further  simplifying  assumptions  are  introduced  in  IID,    in  order 
to  reduce  the  mathematical  difficulties.     In  a  search  for  the  various 
methods  which  would  render  the  general  problem  amenable  to  math- 
ematical analysis,    it  is  found  that  by  defining  a  suitable  mean  temp- 
erature in  the  direction  perpendicular  to  the  direction  of  flow,    and 
by  making  a  further  assumption  of  constant  surface  heat  transfer 
coefficient,    the  useful  formulation  given  in  (II- 10)  and  (II- 11)  is 
obtained. 

Assumption  of  large  velocity  justifies  the  use  of  a  perturbation 
technique,   which  is  presented  in  Section  IV  hereof. 

The  zeroth  order  terms,    which  corresponds  to  a  solution  given 
in  Section  15-3  of  [7]     and  also  in  example  7.  5  of  [2],are  obtained,   and 
first  corrections  to  the  fluid  and  pipe  temperatures  are  obtained 
explicitly.     These  solutions  are  found  not  to  be  valid  at  and  beyond 
the  temperature  front.     Accordingly  a  singular  perturbation  tech- 
nique is  used  to  obtain  an  expression  which  traces  the  behavior  ahead 
of  the  front. 
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In  Section  VI  a  further  physical  assumption,    namely  neglecting 
axial  conduction  in  the  solid,    leads  to  a  problem  for  which  we  are 
able  to  find  the  closed  form  solution  by  using  the  method  of  Laplace 
transformation.     This  affords  an  improvement  over  the  solutions  given 
in  [7]  and  [2],    cited  in  the  preceding  paragraph. 

In  Section  V  the  physical  assumptions  that  are  introduced  consist 
of  considering  an  average  fluid  temperature  and  neglecting  axial  con- 
duction in  the  solid.     This  case  is  solved  using  the  method  of  Laplace 
transformation.     Even  though  the  methods  of  transform  calculus 
reduce  the  partial  differential  equations  to  a  system  which  can  be 
readily  solved  in  the  transform  domain,    sometimes  delicate  math- 
ematical handling  is  required  to  invert  the  solutions  back  to  the  real 
variable  domain.     Such  a  case  is  that  of  equation  (V_27)  and  (V-28), 
where  the  solutions  in  transformed  form  contain  an  infinite  number 
of  singularities   (branh  points).     An  elegant  theorem  for  this  inversion, 
established  by  Professor  J.    E.    Brock,    is  presented  in  Appendix  A. 
Thus  closed  form  solutions  are  obtained.      Furthermore,    use  is  made 
of  the  properties  of  Laplace  transform  to  yield  some  useful  asymptotic 
solutions. 

In  Section  VII  the  view  taken  is  that  of  W.    Munk  [24].     It  was  found 
possible  to  improve  upon  his   solution  by  accounting  for  the  axial 
conduction  in  the  fluid  as  well  as  the  solid  wall.     This  treatment 
assumes,    as  Munk  did,    an  instantaneous  common  constant  uniform 
temperature  at  any  cross   section  of  the  fluid  and  wall.     Although 
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this  result  is  obtained  independently,    it  turns  out  to  be  known;  see 
Appendix  II  of  [13].     This  solution  is  useful  for  small  fluid  velocities; 
Cf.Hiep  [13],    p.   32  and  Preston  [27]. 

We  will  make  no  further  attempt  to  relate  the  problems  considered 
herein  to  practical  engineering  situations.     Such  references  as  [7] 
and  [2]  are  devoted  to  problems  having  similar  relations  to  practical 
engineering  situations.     Their  importance  lies  in  the  fact  that  they  do 
provide  definite  mathematical  solutions  to  well  posed  mathematical 
problems  which  have  relation  to  the  physical  world  of  heat  transfer 
engineering;  it  is  the  obligation  of  the  person  who  employs   such 
results  to  investigate  the  degree  to  which  the  mathematical  conditions 
are  satisfied  in  his  own  physical  application. 

B.      RECOMMENDATIONS 

Most  important  is  treatment  of  the  integral  equation  developed 
in  Section  III.     Remarks  about  this  are  included  in  subsection  IIIF. 

It  would   be  useful  to  develop  asymptotic  formulas  to  employ  in 
place  of  the  exact  inversion  in  cases  where  there  are  severe  numerical 
difficulties  associated  with  the  latter.     The  method  of  steepest  descents 
(e.g.    Cf.    Carrier  Crook  and  Pearson  -  Section  6-3  [8]),  appears 
promising  for  this  purpose. 

There  are  several  mathematical  details  which  should  be  filled 
in  (although  in  treatments  similar  to  our  own  they  are  usually  glossed 
over).      These  all  relate  to  locations  of  the  zeros  and  poles  of 
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transcendental  functions.     In  the  case  of  equations  (V_ 5 9)  and  (III- 155), 
we  have  simply  assumed  that  the  zeros  and  poles  are  all  real;  under 
this  assumption  they  are  not  particularly  difficult  to  locate.     How- 
ever a  proof  that  there  are  no  complex  zeros  or  poles  is  really 
required. 

Inasmuch  as  instantaneously  higher  rates  of  heat  transfer  are  to 
be  expected  at  and  near  a  rapidly  moving  front,    it  might  be  expected 
that  an  analysis  employing  finite  speed  of  heat  propagation  would  be 
enlightening.     Accordingly,    it  is  suggested  that  such  an  analysis  be 
attempted  in  one  of  the  simpler  subcases  dealt  with  herein;   Cf.  ,    for 
examples,    Section  1-6  of  [10]. 

Also  since  numerical  difficulties  are  to  be  expected  in  applying 
the  inversion  theorem  of  Appendix  A  to  the  heat  transfer  problems 
herein,    it  is  suggested  that  actual  numerical  calculations  be  under- 
taken.    A  multiple-precision  facility,    such  as  described  in  [14]  and 
[28]  seems  necessary  for  this  purpose. 

Another  thing  which  comes  to  mind  is  to  relax  the  specialization 
herein  because  of  the  assumption  of  slug  flow,    i.  e.  ,    constant  fluid 
velocity.     It  seems  clear  that  the  assumption  of  velocity  as  a  given 
function  of  the  radial  or  transverse  dimension  would  provide  definite 
differential  systems  for  analysis,    but  that  these  would,    in  general, 
be  much  more  difficult  to  deal  with  than  those  considered  herein. 

The  treatment  given  by  Hiep  [13]   relies  heavily  on  the  accuracy 
and  validity  of  a  numerical  inversion  of  Laplace  transform  due  to 
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Salzer  [30],    [31],    [32].     However  our  experimentation  with  such  methods 
of  inversion  showed  that  reliably  accurate  results  could  be  obtained 
only  in  limited  ranges  and  for  relatively  simple  functions.     This 
persuaded  us  to  look  for  other  analytical  methods  of  treating  the 
problem,    such  as  perturbation  method  treated  in  Section  IV  hereof. 
However,    it  is  probable  that  further  development  of  techniques  of 
numerical  inversion  may  provide  a  useful  and  powerful  tool  for 
analyses   such  as  those  studied  herein. 

We  should  also  mention  our  feeling  that  no  great  difficulty  would 
be  associated  with  inferring  response  to  a  unit  step  change  in  left 
end  temperature  but  retaining  the  other  conditions  specified  in  the 
analysis  of  Sparrow  and  DeFarias  [33],    who  employed  a  sinusoidally 
varying  left  end  boundary  condition. 

Also  although  throughout  we  have  assumed  that  the  conduit  had 
perfect  exterior  insulation  so  that  there  was  no  heat  flux  through  this 
surface,    it  seems  clear  that  other  exterior  boundary  conditions  (such 
as  an  isothermal  condition)  would,    in  most  applications  considered 
herein,    involve  only  changes  in  detail  and  not  essential  changes  in 
procedure. 
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APPENDIX  A 
INVERSION  OF  SOME  LAPLACE  TRANSFORMS 

The  theorem  stated  and  proved  in  this  appendix  provides  the 
Laplace  inversion  of  (V-27),    (V-28),    and  (VI- 14)  of  the  text.     This 
appendix  is  largely  a  paraphrase  of  [6]. 

THEOREM:     If  the  function  to  be  inverted  is  of  the  form 


f(S)  =  _L  ^ 


-l/mcs) 


(i) 


where  m(s)  is  a  meromorphic  function  whose  zeros  s  =  -a^  and  whose 
poles  s  =  -t^  satisfy 


o  <  <L,  <b,  <  Q,<  bi    -  •  •     <%<K  '••  < 


oo 


(2) 


Then,    for      t  >   0,    the  inverse  Laplace    transform    f(t)-=   £.     lf(*)J 


is  given  by 


-it 


(3) 


K«l 


a* 


where 


Haj  =/Fhl 


(4) 


The  interpretation  of  the  sum  in  (3)  is  as  follows.     If  there  is  a 
finite  number  N  of  pairs  a.    ^  b^,    no  special  interpretation  is 
necessary.     If  there  is  a  finite  number  (N-l)  of  pairs  a      <  b     and  an 
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additional  zero      s  =  -  a        aN>bN_i»   an  ^       P°le  s  ~  "  b^  =  -  <*> 

is  to  be  added  so  as  to  make  N  pairs.     If  there  is  an  infinite  number 

of  pairs  a^.    ^  b,  ,    take  N  =  co 

In  proving  the  theorem  stated  above  we  consider  a  sequence     I  R      I 

of  positive  numbers,    R     ^  a    ,    b^  for  any  n  and  any  k,    R  — ►<»  as 

n — >oo,   and  we  consider  a  corresponding  sequence  of  closed  contours 

C     in  the  complex  s-  plane  each  of  which  encloses  the  pole  of  f(s)  at 
n 

s  =  0  but  excludes  all  of  the  branch  points  at  s  =  -a,    and  s  =  -  b^. 
The  contour   C     consists  of  four  subcontours   Cni,    Cn2,    Cn3  and  C    a 
as  follows  (see  Figure  A_l).     Cni  is  a  vertical  line  from  s   =  -  iR   , 

to  s  =        iRnj    indented  as  shown  by  a  small  semicircle  passing  to  the 

*-e      if  w 

right  of  s  =  0.      Cn2  is  the  arc  s   -  Rnf     >   ]l  4  Q  KV       and  C    *  is  the 

similar  arc  s  =  Rn^  ;-Tf^05~-^    •      Finally  Cn3  is  the  contour  which 

passes  to  the  right  froms    =  -Rn,    along  the  negative  real  s-axis  to 

8  =  -  sl\  with  small  semicircular  indentations  into  the  upper  half  s  - 

plane  about  each  zero  (-a^)  and  pole  (_b^)  and  then  passes  to  the  left 

from  s  =-  a,  to  s=  -  Rn  with  small  semicircular  indentations  into  the 

lower  half  s  -  plane,    as  indicated  in  Figure  A_l. 


Figure  A-l:     Contour  C 


n 
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Clearly  by  Cauchy's  theorem 

i    •'«* 


-  \JnUo) 


s    fixp 


-n°) 


for  n  =  1,  2,  3  .  .  .     .     Also,    by  the  complex  inversion  theorem  for 


(5) 


Laplace  transforms 
_      1 


si 


Im-  ^  /*«*>*— *F(t; 


(6) 


fC». 


as  n  — >  00    • 


To  proceed  further  we  wish  to  establish  an  inequality  for  arg 
[m(s)].     To  do  this  we  consider  complex  vectors  emanating  from  the 
points  s  =  "*ai<.  and  s  =  _bjc,    and  we  write 

s  _ 


OK  ' 


(7) 
(8) 
(9) 
(10) 


where 


W 


» 


f»(0*  /i(o)  TT   ^n-rTT  Vk 


(11) 


and 


(■»> 


(n) 


^ = L^  ■  & 


(12) 


The  sums  and  products  written  with  a  superscript  (n)  indicate  that 
only  those  terms  are  to  be  considered  for  which  an  <    Rn  or  b       <     R    , 
as  may  be  appropriate. 


n 
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Now  note  for  any  s  in  the  (closed)  upper  half  plane  we  have 


tk    >  YK  J  K  =  \jl,  •-. 


(13) 


A  geometrical  argument  is  given  in  Figure  A-2.     Now  either 

there  are  N  pairs  a,    ,   b,    each  less  than  R    ,    so  that 

K         k  n 

or  there  is  an  extra,   unpaired  zero,    so  that 

A/-/ 


ft«£(VYj  +  «V>* 


K  =  | 

In  the  first  case  we  can  write 


(15) 


AM 


/»*-  %-  ]T(y*-<U-y, 


(16) 


S=-kk 


S=-a, 


s»r 


Figure  A-2:     Demonstration  that       4*  >  V" 
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whereas  in  the  second  case  we  can  write 

A/ 


A,-  v£(YK-tKti)  $*iS* 


(17) 


k*/ 


Thus  in  either  case 

o  <  p»  <  re  (is) 

Similarly  by  considering  an  s  in  the  lower  half  plane  we  find 

-tf</3„    i°  (19) 

so  that  in  all  cases 

-T7/<  p„  <rr  (20) 

Thus,    if  we  write 


■Lm,  — 


st 


— .      e  f(s)ds 


-*Z 


we  note  that 


<21) 


|™|-    jj;    -p(V^") 


«"  gr   **P  \-fr  e«(Mi) 


(22) 


J     =      Rn 


uniformly  on  Cn2»   by  virtue  of  the  inequality  (20).     Thus,    by  Jordan's 
Lemma  In2   — >0    as  n — +oo.     Similarly  the  integral  on  Cn4  vanishes 


as    n  — ^oo    • 


Next  we  consider  the  integrals  around  semicircles  surrounding 
poles  and  zeros  on  contours  Co        The  result  is  obvious  near  a  zero 


st      -a„t 


s  -  -  ak;    e  &  e     ,   I  ^) 


J-  s 


ince  m  Cs)  ^^  »   and 
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the  length  of  the  arc  goes  to  zero  as  its  radius  does.      Thus  the  integral 

on  such  semicircles  about  a  zero,    s  =  -  a.    vanishes  with  the  radius. 

k 

However,    for  the  integral  on  a  semicircle  about  a  pole,    s   =  -  b^,    we 
must  invoke  (20)  again.     We  have 


<2*p[i/m(7)]   =  e*P  [- 1*„  Cos  (  f~/A     <  i 


(23) 


even  though    M  — *•  oo     as  the  radius  of  the  indentation  goes  to  zero. 
Thus  the  integral  on  such  a  semicircle  about  a  pole,    s  =  -  b^,    vanishes 
with  the  radius. 

We  remark  that  although  we  have  reason  to  believe,    on  physical 
grounds,    the  correctness  of  an  evaluation  by  Kantola  [18],    we  feel 
that  his  argument  intended  to  show  vanishing  of  such  an  integral  on  a 
semicircle  surrounding  a  branch  point  is  deficient.     Our  argument, 
based  on  (20),    is  not  applicable  to  his  case. 

If  now  we  write  (5)  in  the  form 


Li=I     -  I„?   -  I„,  -  I, 


(24) 


Lnl  "  xn  "  xn2  "  ^3  "  xn4 
and  let  n >  oo     ,    and  also  let  the  radii  of  all  semicircular  indenta 


tions  on  Cn3  go  to  zero,    we  find 

F(t)  =  exp  [-  Ji(0)  ]  -  I3 
where 


h  = 


irtL 


jL  £*p  (si-  v/mcsy)ds 


(25) 


(26) 
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and  C     now  is  a  contour  extending  from  s  =  -  00    along  the  upper  edge 
of  a  branch  cut,    to  the  right,    passing  around  s  =  -  a       ,     and  then 
extending  along  the  lower  edge  of  the  branch  cut  to  s   =  -   oo     .     The 
branch  cut  is  the  negative  real  s  -  axis  from  -co    to  -  a,.     There  are 
no  longer  any  semicircular  identations  in  Co. 

To  reduce  (26)  to  a  more  convenient  form,    we  first  introduce  the 
function  U  (x),    for  real  x,   which  denotes  the  excess  of  the  number 
of  zeros   (s  -  -  a^.)  over  the  number  of  poles  (s  =  -  b^)  to  the  right  of 
the  point  s   =  x.     It  is  clear  that   V  -0  in  the  interval  -a^-fi      ^   s     ^    -b^ 
while     V  =  l  in  the  interval    _b,    <  s     <£  _a^.     Writing      /?     now,    rather 
than       fln         ,     since  we  have  taken  n  — >  oo  and  all  zeros  and  poles 
are  included  in  our  accounting,    we  note  that  on  the  upper  edge  of  the 
cut     /3    =         VfT  /2  while  on  the  lower  edge  ft  -  —      yff/2.     We  also  now 


write 


rather  than 


r 


.      Thus 


I    =    -i- 


T 

'-03 


eocp 


St  -p  e*p(L  Vff/i) 


ds  + 


•oo 


At 

-a, 

iiri  i  s 


S-t  -p  e*p(-iV%j]ds 


(27) 


.I       '  •     st 


=  Sri  I  i" e  r*p  ~?  **** (c**/z)  -^p  -?**?  fc*vir/i) 


ds 


-00 


If  t>    =  0,    the  quantity  in  braces  vanishes,    while  if      l)  -   1  it 

becomes 


(28) 
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Thus 


A/ 


P(s) 


ds 


(29) 


If  we  now  change  the  variable  of  integration  from  s  to      1 
obtain 


s,   we 


A/ 


|*(1) 


^ 


(30) 


and,    thus,    finally 


hi 


u 


p(t;  =eocp  Ui(0 


f(n) 


(31) 


The  use  of  the  upper  limit  of  summation  is  consistent  with  the  remarks 
following  the  statement  of  the  theorem,    the  proof  of  which  is  now 
complete. 

In  equation  (VI- 14)  the  solution  obtained  in  Laplace  transform 
form  is  a  special  case  in  which  we  have 


mz(s)  =  <  (s+a,)(s+ai)/(s  +  bi) 


(32) 


where  k     >    0    and       O  <  (X,  <  bi  <  #l  .     The  numbers  a^,    z.^  and 

bj  are  known  explicitly  by  (VI-9,    15,    16).     We  have 


Hio)  =     fo  =  i\Jai&i/bi 


(33: 
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and  the  inversion  is 


poi) 


00 

dn    -   i        <L     A^, 
1        TTJ     1 


f(a) 


(34) 
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APPENDIX  B 
NONEXISTENCE  OF  COMPLEX  ROOTS  OF  EQUATION  (V-36) 

Demonstration  that  such  equations  as  (V-36)  and  (V-38)  have  only- 
real  roots  are  occasionally  called  for  in  order  to  employ  Cauchy's 
Theorem  (and  Jordan's  Lemma)  to  evaluate  the  inverse  Laplace 
transforms  encountered  in  heat  conduction  problems.     Carslaw  and 
Jaeger  [7]  have  some  general  remarks  on  this  matter  (Section  12.  8 
IV  ,    page  324),    making  use  of  Sturmian  theory;  their  reference  is  to 
Chapter  X  of  [15],    Polya  and  Szego  [26]  consider  similar  problems, 
most  of  them  concerned  with  the  zeros  of  polynomials,    but  some.e.  g.  , 
Auf.  172,    of  nature  similar  to  ours.     They  use  various  methods  of 
proof.      Professor  J.    E.    Brock  has  constructed  the  following  elementary 
proof  for  the  case  encountered  in  equation  (V-36).      The  case  encountered 
in  equation  (V_6l)  wherein  the  simple  trignometric  tangent  function  is 
replaced  by  a  generalization  involving  quotients  of  expressions  involv- 
ing J-  and  Y-  Bessel  functions,    is  obviously  much  more  difficult  and 
we  simply  make  no  attempt  to  construct  a  proof  in  this   case,    taking 
refuge  in  the  (weak)  argument  that  the  function  we  encounter  is 
"similar"  to  the  tangent  function. 

Thus,    we  consider 

xtanx  =     — i C  >  o    .,        n>o        o<b<»  (1) 

n-  x2- 
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Let  x    =    a  +  ib,   a,   b  real  and  presume  b  i  0 
We  will  show  that  this  leads  to  a  contradiction 

tan  (a  +  ib)  =  (tana  +  i  tanhb)/  (     1  -  i  tana  tanhb). 
Denoting  the  left  side  of  (1)  by  L,    then 


L  =  x  tan  x  =  (a  +  ib) 
which  may  be  written  as 


tji  Tanafc«fo) 


x  tan  x  =    J '         {_ 

D, 


(2; 


o: 


where 

D     =  cos2a  +  cosh2b. 
Next  considering  the  right  hand  side  of  equation  (1)  and  denoting  it 
by  R  we  obtain 


R  = 


where 


t>2-  (n  +  bl-al)  +4alt>2    I c 


At  this  point  it  is  well  to  observe  that  for  b  /  0,    coshb    ">      1 
While  for  any  real  a,    cosa    >    -   1 
and  thus  we  conclude  that 


D,> 


(4) 


(5) 


(6) 
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Also 


t>z  = 


(M-a  +  b2)-3U0b 


>  o  »>o       j     b  >o  (7) 


2      ? 

for  (n  -  a    +b    )  and  (2ab)  cannot  both  vanish. 

Now  if  we  equate  the  real  and  imaginary  parts  of  (3)  and  (4), 

we  obtain 


afr^lfl-  \*Alflllb  ('*  +  £-al)-tqal £-£(*+&- a") 


D, 


D 


(8) 


by  equating  the  real  parts,    and 


a     b_ 


Dz 


(9) 


by  equating  the  imaginary  parts.      From  (9)  above  we  may  write 


a  A^za    _  _  az   A^ll  b        2  6<Z 


l 


b  D 


(10) 


and  substituting  in  (8)  we  obtain 


1/2 


-{a  +  b  j s  (ir 

Now  considering  the  numerator  on  the  right  of  (11).      Call  it  N,    and 
consider  it  as  a  function  of     0     ;  then 


dS 


a*-(ti  +  bz-al>)  +  lo2  -3ai-n-b2  =  f 


(12) 
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First,    supposing  that    P    >      0,    then  the  minimum  value  of  N  occurs 
for  &  =  0. 

Next  suppose  that      f     <.     0,    then  the  minimum  value  of  N  occurs  for 
6  =  n. 

-  n  (  a  +  t>z)  +  (a  +  fc^  >  a 

Thus  in  any  event       fj  ^   0,    and  since  for  any  real  b  ^  0  we  have 

AUi-lUb      .  (14) 

b  ^ 

thus  the  left  side  of  (11)  is  less  than  zero  whereas  the  right  side  is 

greater  than  zero;  which  is  a  contradiction. 

Thus  we  cannot  have  roots  of  the  form  x  =  a+ib  with  a  and  b  real, 

and  b  ^  0.      Therefore  we  can  only  have  real  roots. 
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